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Introduction 


The  purpose  of  this  report  is  to  present  the  results 
obtained  for  magnetostatic  wave  transducers  under  contract 
number  F19628-80-C-0029  from  the  Air  Force  ESD  RADC/EEA 
Hanscom  AFB,  MA  since  obtaining  the  results  presented  in  the 
Interim  Technical  Report 

In  the  earlier  report  the  analysis,  computer  programs 
and  graphical  results  were  displayed  for  the  flat  field, 
basic  theory  magnetostatic  surface  wave  program,  the  micro¬ 
strip  model  for  surface  waves  and  the  dispersion  relation 
for  generalized  surface  waves.  There  was  also  presented 
the  analysis  for  complex  impedance  for  surface  waves  in 
free  space. 

^  We  here  present  the  analysis  and  computer  program  for 
obtaining  the  dispersion  relation,  group  delay,  magneto¬ 
static  wave  power,  radiation  resistance,  radiation  reactance 
and  insertion  loss  for  forward  volume  waves  and  backward 
volume  waves.  Included  are  the  consideration  of  ground  planes, 
apodization  of  the  fundamental  mode,  multistrip  meander  line 
and  parallel  grating  and  isolated  independent  conductors  for 
a  uniform  current  distribution  without  variable  coupling. 

The  analysis,  computer  program  and  its  use  together  with 
graphical  results  are  presented. 

The  analysis  for  obtaining  the  dispersion  relation  and 
group  delay  for  volume  waves  with  arbitrary  orientation  of 
the  biasing  field  is  presented  here.  The  computer  program 
and  graphical  results  for  cases  of  interest  for  the  generalized 


dispersion  relation  for  volume  waves  including  ground  planes 
are  presented.  Also  included  in  the  report  are  the  results 
of  attempting  to  obtain  nondispersive,  electronically  tunable 
time  delay  elements  by  combining  the  group  delays  for  volume 
waves  of  two  different  biasing  field  magnitudes  in  the  same 
specified  orientation. 

We  also  present  the  analysis  of  incorporating  a  non 
uniform,  hyperbolic  current  distribution  in  the  basic  magne¬ 
tostatic  surface  wave  program.  The  computer  program 
accomplishing  this  task  is  also  presented. 


Basic  Theory  -  Forward  and  Backward  Volume  Waves 

The  basic  theory  leading  to  the  dispersion  relation,  group 
delay,  magnetostatic  wave  power,  radiation  resistance,  radia¬ 
tion  reactance  and  insertion  loss  for  forward  and  backward 
volume  waves  for  a  transducer  consisting  of  a  YIG  region 
sandwiched  between  two  finite  dielectrics  (see  Figure  1)  is 
here  outlined.  We  follow  the  previous  work  on  surface  waves  ^ 
We  start  with  Maxwells  equations 

_  nT  _ 

VxH  =  —  ;  V-B  =  o  (1) 

VxE  =  -  i  V *D  =  o 

and  the  constitutive  relations  in  each  region 

B  =  \iQ  (H+M)  (2) 

D  =  e  I 

where  M  is  non  zero  in  the  YIG  region  only.  We  utilize  the 
gyromagnetic  relation  in  the  YIG  region 

=  -Y  MxH  (3) 

and  retain  first  order  terms  only. 

We  assume  the  time  dependence  of  all  physical  quantities 

to  be  e^wt,  we  also  assume  the  magnetostatic  approximation 

H  =  E  =  E  =  o  (4) 

z  x  y 

ioeE.  =  o 
z 


and  no  variation  of  any  physical  quantity  in  the  z  direction. 


We  have 


3E 


y  “  BX 


9  E. 


=  jo)  B 

x  y 


Also 


B  =  y  H 
x  ox 


B  *  y  H 
y  O  y 


in  the  non  YIG  regions,  while 


I*  11 

B  J  °  \  jy 


y )  \  J  21  u22  /  \y/ 

in  the  YIG  region;  where,  for  forward  and  backward  volume 
waves , 


W12  *  P21  =  ° 


and,  for  forward  volume  waves,  (0=o  in  Figure  1) , 


'll 


(l  +  41Wo\ 


V22  1 


and,  for  backward  volume  waves,  (0=90°,  <{>=0  in  Figure  1) 


y  *  2.8  mhz/oe  ;  4IIM  =  1750  oe  ,  f  = 


where 


We  define 


For  forward  volume  waves 


sin  at k I d  e 


-  k  d 


F(k)  =  s^.n-a  I  ^  I  ^  e -  |-a  cot  a|k|d  [  (coth  | k | t^+1) -e~2 ^ k ^ 

(coth  |k|t^-l)]  +  (a2  coth  |k|  t^-D+e-2  ^k^(a2  coth  |  k  1 1^+ 1)^ 


and,  for  backward  volume  waves 


sin  a  I k I d  e 


-  k  d 


-2  ikU 


sin  a  k  d  e  1  1  (  .  i ,  i  ,  r,  ..  i.  ...  —  2  j  k  £ 

F(k)  =  - — -  <-a  cot  a|k|d  [(coth  [klt^+D-e  1  1 

(coth  Iklt.-l)]  +  (a2-coth  |k|t,)-e  2  ^  k  ^  ( a2+coth  |k|t.)l 

1  1  1)  (2o) 

where  Z,  d,  t^are  the  widths  of  the  three  regions  (see 

Figure  1) . 

If  we  require  that  F(k)  vanishes  we  obtain  the  dispersion 


relation 


tan  alkld 


a  [(coth  fkftj  +  D-e' 


(coth  lkftj^-1)] 


[(a2coth  |k|t^-l)+e  2^k^(a2coth  |kjt^+l)]  1  forward  waves 

|  [  (a2-coth  |k|  t1)-e“2^k^(a2+coth  |  k J  t1>  ]— 1  backward  waves 

(21) 

which  is  solved  numerically,  by  iteration,  to  obtain  k  as  a 
function  of  f. 

There  are  an  infinite  number  of  solution  modes  of  (21) 
corresponding  to  the  multiplicity  of  the  inverse  tangent  func¬ 
tion.  For  each  mode  the  solution  for  positive  k  (s  =  |k|/k  =  1) 
is  the  same  as  for  negative  k  (s  =  |k|/k  =  -1) since  only  the 
magnitude  of  k  appears  in  the  dispersion  relation. 

The  bandwidth  for  which  the  solution  is  obtained  is 
given  by 

Y  Ho  <  f  <  yVho  (Hq  +  4JIMo)  (22) 


6 


We  can  now  find  the  group  delay  from  the  dispersion 


relation  curves,  k  vs  f,  by  numerically  applying 

3o; 

vg  =  3k 

We  can  find  all  the  constants  in  the  solution  and  pro 

f4l 

ceed  to  find  the  field  equations  as  for  surface  waves  . 

r4l 

The  magnetostatic  wave  power  is  then  obtained  from 


f\  r 

2  J  z  \ 


-U+d) 


utilizing  (5) . 

We  then  find 


P  =  a  po  G  A 

4|k|2 

where  k  is  obtained  from  the  dispersion  relation,  and 

_  - 1  k  I  d  = . 

G  =  e  1  1  J  (k) 


3k,F(k,)  k’=k 


For  independent  conductors 


[  J  (1c)  J  =  |  2  sinc  w~  n1 


where  a.^,  p^  are  the  strip  lengths,  strip  widths  and 

center  to  center  spacings  for  the  N  strips  and  n=l  for  a 
parallel  grating  while  n  =  -1  for  a  meander  line  and 


sine  x  = 


sin  Ilx 


For  uniform  N  strips  we  can  write 


N  jkpN 

|S(Ml  =VVsinc  jfrl  I1'"  fkp 

1-n  eJ 


For  the  quantity  A  in  (25)  we  have,  for  forward  volume 
waves , 

A  =  2e-2lki*(sinh  2 1  k  i  i-2  |  k \l)  +  <a2+l)  d  |  k  |  -  Sin  2a!k!d+cos  2°t,k|d’1j 


t  [sinh  2lkl  V2lklMraVl  .  jo^ll  cos  2a | k | d  -  a  sin  2a|k|d] 
2  sinh  | k i t ^ 


2e-2  I k I  ^  |(a2-l)d|k|- 


(a2+l) 

2a 


sin 


2a|k|d  +  Uinh  2Wy2|K|y 


2  sinh  Ik)^ 


t2?=I  _  cos 

L  2  2 

2 


2a  |  k  |  d]| 


(30) 


+  e 


— 4  |  Jc  |  £  j(a2+l)d|kj~  sin  2a  |  k|  d  -  (cos  2a|k|d-l) 


fsinh  2 1  k  1 1,- 


tj-2 1 k) t A  ca 
2  lk|t~  / 


L±I  _  -L  cos  2a|k|d+2  sin  2a|k|d]| 


2  sinh  ]k|t^ 
and,  for  backward  volume  waves, 

A  w  2e-2,kU(sinh  2|k|l-2|k|n-(u2tl)|k|d-^ii  sin  2c.|k|d+cos  2»|k|d 


/sinh  2|k|tr2|k]tl\  /2+1  ^  ^  ^  2a|k|a  +  sin_2|ikid\ 

+  ^  2  sinh2  |k|t ~/W  ^  ' 

+  2e'2|k|l[(cl2-l)d|k|+  sin  2a|k|d  -  /81nh  2|k|t^  2|kl^ 

\  2  sinlr  Jk|  tx  / 


(  ill  +  2-ii  cos  2a  |  k  |  d  )  3 
2aZ  2aZ 


(31) 


+  e_4lki  a^l-(a2+l)d|k|-  sin  2a|k|d  -  cos  2a|k|d 
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fsinh  2  |k  1 1^2  |k(  tjj 
2  sinh  Ik  it. 


|0.2+1  +  JLzl  cos2a|k|d  -  — n  2alkld' 


2  a 


2  a 


a 


Proceeding  as  for  surface  waves^1^,  the  radiation  resis¬ 
tance  is 

4  |  P  | 


R  (1-n)  +  (l+n)N2 

for  one  wave. 

The  total  radiation  resistance  is 

R  =  2R 
m 

and  the  radiation  reactance  is  obtained  from 


(32) 


(33) 


xm(£)  -  I  / 

—  oo 


R  (f')df' 
m 

— — 


(34) 


by  numerical  integration. 

The  insertion  loss  is  then  given  by 


IL  =  20  log 


10 


4  R  Rg  _ 

(R  +R  +R_)2+(X  +Xt)2 
g  m  L  m  L 


76.4  10* 
9uj/3k 


AH  Ar 


(35) 


where  R  is  source  resistance,  R_  is  conduction  loss  and  X_ 
g  Jj  l 

is  a  series  matching  impedance.  AH  is  linewidth  representing 
material  loss  and  Ar  is  the  propagation  distance. 

We  thus  have  relations  for  all  the  quantities  of  in¬ 
terest  for  both  forward  volume  waves  and  backward  volume  waves. 

In  figures  8-17  are  presented  graphical  displays  of  the 
dispersion  relation,  group  delay,  radiation  resistance, 
radiation  reactance  and  insertion  loss  for  forward  volume 
waves  and  backward  volume  waves. 


9 


Dispersion  Relation  -  Generalized  Volume  Waves 

In  this  section  we  obtain  the  dispersion  relation  for 
volume  waves  including  ground  planes  when  the  biasing  field 
can  be  in  an  arbitrary  direction  (see  Figure  1) .  We  have 
the  permeability  tensor  components  appearing  in  (7)  as^^'^5^ 


Mii  =  1  +  y2H  (4IIM  )  (sin2  0  sin2  <p  +  cos2  0) 

11  o  o  - 5 — 2 - 2 - - 

y  H  -  f 
o 


U,_  =  1  +  y2  H  (4IIM  )  sin2  0 
li  o  o 


"2 — 2 - T 

Y  H  -  f 
o 


- j Ml2  =  jY  (4IIMo)  sin  0  (f  sin  4>+jY  Hq  cos  <J>  cos  0) 


(3fc 


~2~J2 - IT 

Y  H  -  f 
o 


jM21  =  -jY  (4IIM<_i)  sin  0  (f  sin  <|>-jY  Hrt  cos  <|>  cos  0) 


— 5 — 5 - T 

Y  n  -  f 
o 


where 


U> 


Y  =  2.8  mhz/oe  ,  4IIMo  =  1750  oe  ,  f  = 


(3 


The  solution  is  in  the  form  of  a  potential  function 
such  as 

^  =  ejUt-kx),Re|k|y  +  Be-|k|y)  (3; 

in  the  non  YIG  regions  while,  in  the  YIG  region,  we  have  the 
form  for  the  solution  as 


_  e j  (wt-kx)  cqs  c | k | y  +  B  sin  c | k  jy) e"*^kby  (3' 


:2  -  -(u21-u12)  -  4l,ll  v22 


where 


and,  for  volume  waves, 

c2  >  o.  (41 

Here  we  note  the  additional  propagation  term  containing 
b  in  (39)  which  vanishes  for  standard  forward  volume  waves 
(0=o)  and  backward  volume  waves  (9=90°,  <{>=o)  . 

By  attempting  to  satisfy  the  continuity  conditions  and 
boundary  conditions  we  obtain  the  dispersion  relation 

[ (coth  | k | t^+1) -e~2 ^ k ^ (coth  |k| t^-1) ]cy22  cot  c|k|d  = 
coth  |  k  1 1^  [c2  y22  —  ( j  1  +  bV>22)2^ 

-1  +  j  (coth  | k 1 1^-1)  (jP2i  +  by22)  <42 

+e“2k<l|coth  |k|t1  [c2  u22  -  (ju2i  +  by22)2-^ 

+  1  -  j  -Lj^-  (coth  | k 1 1^+1)  (jy2i  +  by22^ 

There  are  an  infinite  number  of  solution  modes  for  this  equa¬ 
tion  corresponding  to  the  multiplicity  of  the  inverse  tangent 
function.  For  each  mode  there  is  one  solution  for  positive 
k(|k|/k  =  s  =  1)  and  another  solution  for  negative  k(-I-£-L  -  s 

This  dispersion  relation  reduces  to  the  ones  obtained 
for  forward  volume  waves  and  backward  volumes  by  specializing 
the  6  and  <{>  values. 

The  bandwidth  for  the  existence  of  volume  waves  is 

[2] 


known  to  be 


The  dispersion  relation  curves  obtained  in  the  plane 
defined  by  4>=o  consist  of  a  backward  volume  wave  curve  with 
a  forward  volume  wave  curve  as  shown  in  Figures  5  and  6. 

The  frequency  at  which  the  changeover  from  backward  volume 

r  2 1 

wave  to  forward  volume  wave  occurs  is  given  by 

fc  =  yYh0(Ho  +  4ITMo  sin2  0)  (45) 

This  critical  frequency  is  at  the  midpoint  of  the  frequency 
spectrum  when  0  is  approximately  43°  as  is  shown  in  Figure  6. 

It  then  appeared  feasible  that,  if  two  volume  waves 
with  different  biasing  fields  were  joined  at  this  angle  0, 
we  may  be  able  to  obtain  regions  of  constant  time  delay. 

Figure  7  shows  the  results  of  such  attempts,  one  curve  is 
the  joining  of  volume  waves  with  biasing  fields  of  2500  and 
2125  (oe)  and  the  other  curve  is  the  joining  of  volume  waves 
with  biasing  fields  of  2500  and  2200  (oe) ,  all  at  the  same 
position  <j>=o°,  0=43.19°. 

The  results  obtained  can  be  compared  with  those  ob¬ 
tained  from  the  joining  of  a  surface  wave  with  a  backward 
[3] 


volume  wave 


Hyperbolic  Current  Distribution 

To  be  considered  here  is  an  alternative  to  the  flat 
field  current  distribution  J{x)  which  will  now  lead  to  a 
replacement  for  the  transform  function  J(k)  of  equations 
(27)  or  (29). 

With  l a,  p  being  the  strip  length,  strip  width 
and  center  to  center  spacing  for  the  N  transducer  strips, 
the  current  distribution  now  under  consideration  is 


i=l,2,...  ,N 


.  .  I  cosh  -x~.(^~.1l£  i_i 
J(x)  =  _ o  6  n 

a  "  sinh  a/2<S 

a/25  -|+(i-l)p<x<|+(i-l)p 


(46) 


(see  Figure  2) . 
where 


°1 


(47) 


for  some  given  6.^  ;  j  being V-l  . 


With  1Q=1»  the  Fourier  transform  gives 


|  J(k)  |  = 


■V£1 

/l-n"  ejkpN\ 

sinh  a/26 
a/26  “ 

\l-n  e^P  J 

[Aje 


a/26+A2e-a/26-j(A3ea/26+A4t 


a2  (k2  +  j  ^) 


where 


A^  =  dj^  (cos  c1  +  cos  f^)  +  f^  sin  c^  +  c^  sin  f^ 


(48) 


A2  s-d^^  (cos  c^  +  cos  f ^)  +  f^  sin  c^  +  c^  sin  f^ 


(49) 


A^  =  d^  (sin  -  sin  c^)  +  f^  cos  c^,  -  c^  cos  f^ 


d^  (sin  f^  -  sin  c^)  -  f^  cos  +  c^  cos  f^ 
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COMPUTER  PROGRAMS 


A.  Basic  Theory  -  Forward  and  Backward  Volume  Waves 
There  is  a  computer  program*  operational  on  the  CDC 
6600  at  Hanscom  AFB,  MA  which  incorporates  the  results  of 
the  basic  theory  for  both  forward  volume  waves  and  backward 
volume  waves.  The  program  produces  plots  on  the  Calcomp 
plotter  at  Hanscom  AFB  depicting  the  various  physical  quan¬ 
tities  as  functions  of  frequency.  Plots  as  well  as  print¬ 
outs  are  obtained  for  wave  number,  group  delay,  radiation 
resistance,  radiation  reactance  and  insertion  loss. 

Flexibility  is  designed  into  the  program  so  that  one 
can  choose  the  solution  mode  for  (21)  and  whether  we  have 
uniform  conducting  strips  or  apodization  in  strip  length, 
strip  width  and/or  center  to  center  spacing.  Additionally, 
the  program  provides  the  relevant  frequency  range  from  (22) . 

Now  follows  a  detailed  description  of  the  input  cards 
with  details  on  the  use  of  the  above  described  features. 
Columns  1-72  may  be  used  on  the  first  5  cards  and  columns 
1-70  may  be  used  on  the  last  three  cards. 

Card  1-Hq,  t^,  d,  JL,  mode  number,  option 

These  six  quantities  are  here  supplied  separated  by 
commas.  Lengths  are  in  meters.  Mode  number  is  o  for  the 
fundamental  mode  of  the  dispersion  relation  (21).  Option 
is  o  for  forward  volume  waves  while  option  is  1  for  backward 
volume  waves. 

♦FORTRAN 
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Card  2  -  AH,  Ar,  N,  n,  RL 

Here  are  entered  these  five  quantities  as  defined 
earlier. 

Card  3  -  first  A£.^,  option 
Card  4  -  first  a,  Aa,  a  option 
Card  5  -  first  p,  Ap,  p  option 

In  cards  3,  4  and  5  three  items  are  entered  for  each  of  the 
quantities  a,  p.  For  each  card;  if  the  option  (third 

item)  is  o  then  the  dimension  is  entered  for  the  first  strip 
(first  item)  and  an  increment  is  entered  for  the  remaining 
strips  (second  item),  but  if  the  option  (third  item)  is  1 
then  the  dimension  is  entered  for  the  first  strip  (first 
item)  and  the  increment  (second  item)  increments  the  next 
strips  and  is  then  a  decrement  for  the  last  strips. 
The  number  of  strips,  N,  must  be  odd  if  the  option  is  1. 

When  there  is  no  apodization  the  increment  (second  item)  is 
o  as  is  the  option  (third  item) . 

Card  6  -  heading  for  plots 
Card  7-  heading  for  plots 
Card  8-  heading  for  plots 

These  three  lines  will  appear  as  headings,  in  the  same 
order,  on  several  plot  frames. 
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F3CG'*4M  VCl  IN  C  IlSPU ,  CUTFUT  ) 

OPN  =  f!  *‘--F  C  <.WA RL  VOLUME  WAVlS:5"  OPN  :  1  «--GACKMRC  VOLUME  WAVES 
C1M.KETCM  P*CGIC  (3),F  l  12C1>  ,FP  (12011  ,FM(  12C1)  ,CAP  ( 32GD  , 

XL  rt*  (IE  Ul)  ,  VOPC  12  G1  >,  VIM  (1201 ) 
c  IM? NCI  ON  Hi  AJ (?) t  HfcAUl (7)  ,HEAC2  <  7) 

C iWf KilCM  F  P ( 1 2u 1 )  *PM ( 1201)  ,  w>F  (12  u  1  >  ,  FM  <  3  2  01  >  »RT  (  1 2  Cl  > , 

X  P  X  ( 1 2 1 1 )  »  j  l'  '  P  ( 1 2  3 1 )  ,SE^M(12-1» 


CCMWCN 

P  -  t.  0  “ 

n  .  a  i;  * 

'  A  i  ’ 
i  AC  ♦ 
**  ;  t  .1  * 
t  h  /« i 
A  fc  A  i.  1 
e  -  A  1 


IN  ll.hl.Tl  ,S  ,C  ,UY  Y  ,PI,  ENM  ,  FCA  ,OPN  ,E  T  i  ,EN,tll  (30  ),FE(33  ),  AA  (  2C> 
“.H,Tl,0,FL,EKM,OPN 
*»□;. L  P,01ST,tN,t  I  A.RL 
♦,El  OtGNttLCEL.ELCFT 

♦  ,APF.GiN,  A ucl  ,A(1P  T 

*  ,  P  ai'  G  T  N  *  PC  EL  , POP  T 
1  i)  u  «  H  c  \  Li 


1  U  u  i  "  i  -u; 

lOO.Hf ACi 
131,HF432 


LC  Ml  1=1, N 

Ell  (I)  =  EL3fGNt  (1-1  >*ElQcL 
A  a ( 1 ) =ASCG INM  1-1)  *AOt  L 
Fl (I)=PPFGIN+(i-l)*PUtL 
N‘l  =  (N  +  l>/2 

IP  (CLnPl  ,,:o.  t.)  GO  TO  hZ 

rr  1=  Nr'l  l  ,  N 

r  I :  (I)  =  r.  LI  (NEl  )-  (I-N£l  )*ELOEL 
1 F  (  A  OPT  ._fl.  j .  )  GO  TO  «*<♦ 

CC  -S  I  =  NcL,N 

AfU)  =  Aa<NLI)-(1-NEL>  «acf.l 
IF  (FCPi  .  E 0.  Cl.  )  GO  10  Afa 
Z C  mV  1=MFL,N 

P'l(  I)=P£  (Nt  D  -  (  T-FFL)  ♦  PC£l 

CCN t TKUE 
F ' I  NT  no 

P  'IN’Tel  ,  H ,  TlO  , f  L,ENM 

F'lNT  *,**  TcLTA  h  =  ”  »OELH,”  DISTANCE  s’SCIST,”  CFTION  =  ",OPN 
P-TM  *,  "  P  LCSS  IS  ’■  ,  ;L,M  N  =  ETA  =  *‘,FTA 

F-iM  f.l*(I*tL3  (I)  ,1  =  1  ,N) 

F-]M  12,  (  1,  AA  ( I)  ,  ]=t,  N) 

P  »  N  '  A3,  (  I, PE (  1)  ,1  =  1 , N) 
u  >•  TNT  c- L 

l;TA  FFCCl  L/tMPEIN?r.RG  ,AH213m,cUVOLUWF/ 
cm  FfcTI'Jo  (PvCGIC,2JJ  .,32..:.  ) 

«=  1=  F.iAlL'.JcE-. 

NC=CFP 


SCN=  (-1  .  >*■  *  NO 


P  h  T  -  .  . 

iCFN  .  ; .  f. ,  J.  )  T  P  ST  =  ')  . 

IF  (CFN  .CC.  1.)  TPST=90. 

^■^l-THCT-»K/icO. 

t  ll  =  v. 

p  r-  -  •  :  -  . 

Lr  =  t..*PI‘»l 
1 NP=EL 
■ -l  =  T  1 

1  '  •  I  HP 
■  N  -  1  1  i  fc  . 
r  .■  Ms?  ,  i- 

r  .  r--(  .3  '•flGFT  (H'  <  H*l  vr?‘j  .*  (SIN  (T  PET  I  )  **  t*  (SIN  (FH  )  )  *'•2)  ) 
Fr  1-2  .  i  GO  NT  (H1  (  P+  1V5  L  •  )) 

F  .t  i 

f  -i  G  =  A  [NT  (  Pi  C)  +2C. 

*  F= INF (FHI ) -INT (FL  C) -  3  ft 

Nf=NF/2 

uC  :  i  -  i  ,  Nf* 

f  U  >  =  f  *t.G+  ( T-i  ) 4  F  0  EL 

T!  lES^-r^K 
f  -  1 

l 

1 

1  f  =  f  ( «  ) 

1  C  (GAM*  h)  ♦*  c-u  F  ? 

IVY  =  IL/If.  *  (SIN  (T  PET) )  ♦  *2*1  , 

fc>X  =  iL/i  Z*  ( (Si  N  (  TPF  I)  *  SIN  (PHI  >>--2MCCSnPET>>**2>41. 
T'M=  -  ,  •  *  L  X  x*i|  i  Y 

IF  ((EM  .lr.  o.)  P  F  IN  I  NEGATIVE  SOLAFe  fiCCT" 

lr  (TEN  ,L  E.  t,  .  )  GO  TU  35 
fc  .  OFT  ( i  EM  )*  ,5/UfY 
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R, produced  From 
best  available  copy. 


If  ( L  .  E  0 .  2)  CO  1C  ? 

HcMo  3 
E  =  -i. 

r  cm  iKur: 

if  <J  .  C-V.  1*  GC  10  b  6 
t  ,;0  =  ATAN  (2  .-r  UYY/  ((C^UYY)  **  ' 

If  1C  2  0  .11.  .)  c  AO=l  AO  +  A0ii  (P  I/C/0) 

C;-0  =  C20  +  i!;ii.(  :,iM*PI/C/L  ) 

I  CM  INLu 

if  (I  .i!C.  D  GC  TC  61 
TF  ( J  .  .  C .  i  )  C-C  TO  51 
I  c  <  L  . EQ .  1)  C AC=CAP( 1-1) 

If(L  .EG,  Z)  0  AO  =  c  AM<  j-1) 
f  :  1 

CCNTINUcI 
C  0  J.  =  F 1  f  C  /.  0  ) 

PM M  « ,CAC,CAi ,FC« 

IF  U-S  (  (CAl-t.  AC )  /  C  AO  )  ,LT 
f /-C*CAl 
F  -  M  ♦  1 

X  F  .Gi.15)  GO  TO  13 
GC  TC  5 

T.  r  ft  « 1 0  •  <i)  GO  IC  2j 

C AsCAl 

If-  ( A  3S  <  FC  A  )  .GT.  1.  )  GO 

IF  (Ca  .  IT  .  X  GC  TO  35 

FF  ( I  )  =  LF 
C  TP  (I  >=Ca 
:  =  ]  ♦: 

P-IM  +  « **C  »*=  ",CA,*f=  **  »EF 
1-2 

ro  ic  3n 
ca=ca: 

If  ( i  F  S  (  F  C A  )  .GT.  1.)  GO 

IF  (CA  . LI  .  0.  >  GC  10  35 

FF <  J)=£F 
t^M  t  ))=CA 
J-J  +  l 

0  ■  I N T  ♦,”CA=  •■tCA,"Fa  ",EF 
IF  (J  .10.  I)  GO  TO  15 


■l.n/C/o 


w  OP  1  ) 


TC  1C 


GO  TO 
35 


GO  TO 
35 


IF  ( J  . r- 
I=J-.l 
K=K-1 
J=  J-1 
k  =  K  +  1 
IF  ( K  .1 
P  :  I  K T  t.< 
11=1-1 
Jl* J-1 
GC  IC  2  s 
P.-IM  *  , 


•  It# 

tj 


GO  TO  15 
GO  TC  31 


GC  TO  5  0 


P-'IM  *  ."lit.  RATION  OOtS  NOT  CONVERGE. F  =  -  ,  £  F  ,  *'  S 
If  (l  .cO.  2)  GO  1C  1, 

l  =  2 


Gf  TO  2 
GCNIINUr 

P  "  I  M  1 3  »  (  F  °  ( L  l )  ,CAP(LL).LL=1,11,5) 
P-IM  c-«  (  FN(..L  )  ,CAHUL)*tL-l«  J1  *  5) 
P  ■,  I  M  o  } 

CC  21  I  a  1 1 I L 

IF  (I  .LO.  1)  VGF(  I)*?  ./PI+(CAPf2>-l 


IF  (I  .t‘0.  1)  VC«F(  I)*?  ./PI* (CAP(2 > -CAF  (1  >)/FlEl 
IF  (I  •  t.  C.  ID  VGP  (1)*5./?I*  (CAP(Il)-CAF  Ul-1)) 
IF  (I  .NL.  1  .ANT.  I  .NE.  ID  VGP (I)  = 
r5./fIMCAPim)-CAFtI-l))/Fucl*,5 

C  r ml t i  1 1  r 


1  ))/FCEL 


<5»/fl*(CAP(I+i>-CAF;i-l))/Fucl4.5 
CONTINUE 

PFIM  3/,  (FP(T)  ,  VGFII)  ,1  =  1,11,5) 

DC  22  J=1,J1 

IF  (J  .CG.  1)  VGH(J)  =  ?  ./P1*(CAM2)-CAM  1M/FCEL 
If  U  .EG.  J II  VGM  <  J>  =5./PIMCAH(  JD-CAF  (  Jl-m/FCEL 
If  <J  .Nt.  L  *  ANC  •  J  .Nt.  Jl)  VGM  ( J )  = 

*  5.  /  PI*  (CA  M(  Jti  )  -C  AM  J  -1) )  /FJEL*.  5 
CONTINUE 

P  'IM  35,  (  FT  ( J  )  ,  VGF  (  n  ,  J  =  1 ,  Jl,  5) 

s  *  1  • 

■  ’C  c;  I  =  l ,  I  x 
>  I  t  FP  (I ) 


l  p  -  r.  a  t- 1 1 » 

T3=  <GAM4H)  **■ 
GYY  =  K/  TP*  (S. 
t..XX  =  TU/T  <  ( 
T  c  Mr 

C  =  SC»T  (TEM  M 
02. 4F  I4EF41 
IF  (  C  F  N  .  E  Q . 
TF  (OPN  .EG. 
r 1  =  SGIs4 
X  IC~*2*  J.  . ) 
(C  OTH  (C  *-‘ 

<  . •  (C**r- 1. 
*-  r "Z (2.-c o 
IF  (:'  .  gi  .  :■• 

i  -  ? .  4 

<  •:  xp(-r. 

XS in (2. *C*C A4 
V  (.  r4  <C“4  2-1 
P,  =  i:XF(— ..  4C 

< 

X  L  C  F 
v  -TON  ( 
rc  :c  5f> 


2~£F442 

IN  (T  FfcT )  M42*l  . 

SI  N  (  1 H  E  I  »  4  GIN  (PHI  )  i**  2  4-  (CCS  (1  F  E  T  )  )♦♦  2  >  ♦  1 . 
-  +  .  4  U  X  X*U  V Y 
. ;/ UYY 

*Jl.  )  CN=1. 

i  .)  n n=i . /c+*z 


(C4*2-l.  )  /<!.  /L*aTN(2.  «(+CA«C)  +CCSU  .4C4CA4D) 
r  u)  -  CA'cL  MC5CH  (CA  *l.L  )  >**21  *  1 .  5 4  < C4 4 2  4 1  .  > -SGN*  ( 
>m;oS<2.4g4CA40)*..  4SIN(2.4C4CA40>))40N 

TO 

37.)  GC  TO  be 

*r  A4  T  1  >4  C-2.*CA*Tl  +  (C**Z-l.  )4E4CA-SGN4(C442+l.)/t./C4 
0)  4SGNMC0TH(CA*EL  >-CA4F.l4(CSCMCA4EL  » )442>  4 
.)  -SIGN4.!,4  (C44  2  +  i.  )4  C  OS  (2ivC4CA4C  ) ) 4  C  M 
A  ~  7  i )  4  (SGN  4 

(C442*i  .  )*Q*CA-.54  (C‘4  2-1.  >/C4SlK  (2.  4C4GA40»  - 
C>  ♦  (CCTH(LA4EL)-Ca<’EL4(CSCH(CA4EL))44S)4(.54(C44241  .) 
. G4  <C442-i. )  4C0S  (2,  4G4CA4C)-C4S  IN  (2  .4C4GA40)  )  )4ON) 


ii  INGE 
•  -  £  E  ( (  A  ) 

T):  U*iJiV  G  EG  A 

I)  -  fi  If.  (  HP  (I  I  )  /  (  (X. 

y  J  MJE 


4424  (Pl*f2*P3  +  FA)  /  a  ./CAJ!42 
£IA)+(i.*ETA)4Eh44e»4A. 


X  (. 

P  ..  : 
X 

x  i  r.  s 

x 

r,c 

p  •  = 


cr  w-i • Ji 

KM  J ) 

L  A  P  ( J  ) 

(Gf  M*H>  4*2-FF  44  2 
- 1  G  /  1  "!*  <  3  IN  <  T  P  F  1 )»  4  4  2  +1  . 

-  Tl.'/T!  *  (  (  GI  N  (  Tl-E1  )  4SIN(  Phi  )  M42MCCS.  (TPET)M42M1. 

.  4  LX  X*IJ  Y  Y 
O T  (TEH)*  .G/UYY 
1 .  4P  I4i  F  41 .  '  t. 

(cfn  ..:c.  ,.j  on~i • 

(CFM  ...  f.  1  .  )  C  N- 1 .  /C4  4  2 

FCN4 

(C4  1  .)  4  -'Ca-  (G44  z-l.  )  /  2.  /C*  SIN  (c.  4  (,4CA4C>  +CCS  (2.4C4CA4D) 

U  CTH  <u  a*  El  >  -C  A  *EL  4  (CSCH  ( C  A4.rL  )  )  4  4  l )  *  ( .  5 4  (  C4 4 2  4 1 . )  -  SGN  4  ( 
Mr»<;'-l.)‘COSC,.*C4CA4OMC*SIN(2.4C4eA4C)>)4CI^ 

•'?  (  2 , 4 G  A '  T  1 ) 
t-  .01.  327.)  GC  TC 

*  4 

‘  XP  (-2.  Tl)4(“2.4CA4T  1*(C44«”1.  )4C4CA-SGN4(C442+l.)/2./C4 

I  ( F  «  4  C  *  C  HJ'C>+bGN4lCUTH(CA*rlL)-CA4EL4(CSlF(CA4tL  ) ) 4  4  c )  4 
E  4  (C,42-l  ,)  ”SGN4,54  (C442+1.  )4COS(2»4C4C£4C>)4OM 

i  XP(-£  .  'C  A  1  1  )4  (SON4 

«C442U  .  >4l  VGA-.!M  U.  44  2-1.  )/C4SlM2.4C4CA4:i>  - 
(2.**-1  7  A "  .3 )  ♦  (COTH(GA4cL)”CA4EL4  (CSCF  ( l  A  4 1 L  ))44£>4(.54(C44?41.) 
CN-M  ,54(C44c-1.)  4GCj  (2 . 4C4  C  f4  C  )-C4SIN  (2  ,4C4CA4J  )  )  MON) 

f  0  *•  s 


F  -=  .  . 

C  N  m  U  r 
i;  r:  Arci  ( g a  ) 

»"<J)=  0*00"  CEC.a  4424  (  P1+ F2+P3+FA) /A./GA442 

:  P  ( J  )  -  A  ^  f  (FM  (J)  )  /  (  (l.-fcTA)  ♦  (1,  fET  A)4EN442)  4A. 

f  C  N  i  I  MJ  r 

r  I  *T  7l  .  (  FP  (I  )  ,PF  (I)  ,1*1,1;,...) 
f’-ur  72,  (FM(J)  (J)  ,  l  j ) 

IF  (11  .ME.  J 1 )  C-0  TO  <sO 

:n  Pa  I-i  ,11 
K  )([)-■•"<  i  )  ♦  CM  (  I  ) 

F  •  I  M  rj  ti 

F-IM  75,(FP([),itl  (I),I=1,I1,1J) 
r FI -FT  Cl) 

FFL  =F F  (II  ) 

U’ U  FT^AN  (  <T, px  ,  1  i,fp i,fpl ) 

F'lM  c  0 

p*  11-1 


i; 


P  '  I  NT  77,  (F‘M1  J  ,PX  (1)  ,  1=1, M  ,20) 

C r  «„*-  i  =  ;  ,  m 
JH=.?.--pj  +  F"(I)*ELL 
X  L  -  / 1.  *  5  .  L‘  > 

L  =  I 

IF  (I  .it.  1)  1=2 

1 J  F  p  ( I )  =  2  J .  *  «  L  OG  1  C  <  ( 9  .*KP  ( 1  >  ✓  RG )  / 

X  (  (i.  +  (FT  (DfPU/fift)  ♦*<: 

X  f  (  t  PX  (I)  +XL  )/  ?G)*  '  l)  ) 

S;  'P(I)=S.tP^C  )-7t.u*iJELH+0IST/ 

X  <?.  .c  t  (F  P  (I*  1)-F  P  (L-l  ) )  ✓  (CAP  (1*11  -CAP  (L  -1 )  )  ) 

? .  ^  M  < I ) =  2  0. "AirGl &( ( 9 .*fiM ( I )  /  KG  I  / 

X  ((!.  +  ( :’T  (I  » +•  P* L  )/9G>-  *  2 

X  +  ((PX(I)+XL  4  c)  ) 

S~Mn)=S_‘-M(F)-7t.4*lit.lH«'DIST/ 

X  (2.  '  PI-  1...  fc  (Fh(I+l)-hK(L-l))/(f,AH(In»-(AMl*l))  ) 
CCNTINUt 
P  ■*  I  N  T  oJ 

PrTKT  77  ,  (  F°  (j.  )  , SE  FP(  1)  » 1=  1  ♦  H  ,!«*  ) 

FAINT  cu 

FSIM  7u  ,  (  FP  (I  )  ,  S£  KM(  i  )  ,I=1,M  .10) 

F  =.;C0=.  il^F  (EG 

yviN  =  £TNT  (  F  i  E  G  0  »  ♦  1  CL. 

XX  =  MM  (  .j  1-FHI)  -  AlNf  (  .ul*FL0>  41. 
c  x  -  2 1  . 

YMI N  =  c . 

C  Y  =  l.Lfc 
CC  12  1=1, II 

IF  <CAP(1»  .  G  i  .  1 ,  F  h )  CAP  (  I )  =i  ,E  o 

i/CP  (1  )  =  4"b  (i/G°  <1  )  > 

TF  <VCP(1)  .  G (  .  UCu.l  VGP  ( 1)  =  l:ju  J  . 


TF 

(VCP  ( 1  ) 

.  G  ( 

LC 

1*.  J=  1 , 

j  1 

IF 

(CAil  ( J ) 

.  Gl 

V#  C  M  ( J  )  =  A  c  S 
IF-  (VCM(J) 
7  =  11. 

X  =  A  *  , 

v  Y  =  j  Cl  . 

Xx  =  _P  . 

COIL  FlC’  ( 


.  G  |.  l.Eb)  CAP  (  J  )  =  1  .Fb 
('JGM  <  J  n 

.  G r  .  lfJCU.)  VGM(J)=  10-10 


CALL  FLG1 
CALL  CYM" 
CALL  CYMP 


CALL  Ftm  (1.5,  c  .,-3) 

L  j  L  L  SY”  40L  (  .  -  t%«6  *.1*HEAD*9*7  3) 

CALL  SYMBOL  (  .  j  ,V.c  ,.l  ,  HE  ADI  ,  0  ,  7u ) 

CALI  EYM-10L  (  .7  ,i  .A  ,.l  ,HEA02 ,0  ,  70) 

CALL  AXK.(v.,  ..2EPWAVE  NUM6E <  t  4- )  (l/MtTEfc)  »26,>*‘=0.,YMIN,CY»YY) 
CALL  AXI?(C  IS  FFR  t  GUtNC  Y  (MHZ), -1**2  ,  0  .  ,X  M I  N  ,C>  ,  t  Y ) 

CALL  LINE  IFF,  CAM, 

Jl,i,U,l,XMlN,CX,YMIK,bY,.C«> 

CALL  F  IGl  (.fc.  ,  l  •  .  “  3 ) 

CALL  CYM"OL  (  .b  ,9  .o  ,.l ,  HEAD  ,  0,  7  C> 

CALL  CYMPO  t ( .I ,9.6  ,.l , HEA01  ♦  J, 7j ) 

C  ALL  SYMTOK ,v ,9.4  ,.l , HEA02.Q,  71) 

CALL  A  X I L  <  C  •  - .  »  2  0  FW  A  V  E  NUMPt*  K  ( -  )  ( 1/ M  c  T  6  K  )  ,2  o  ,  >  ,  9  C  . ,  YMI N,  DY  ,  Y  Y  ) 

FAIL  AXIM  C.  .  ,  lEFF^tQLENCY  (MHZ), -1C, 2  ,  0  .  ,XM1  A  ,LX  ,  Y  Y) 

C-L  L  L  INF  (  f-F.CAF, 

J1  .1,  a  ,1,XMIN,  CX,  YMIK  ,0Y  ,  .Ci8) 

r  y  =  j  j  . 

LALLFLCI  (  1  'i . ,  0  .  ,  -  3 ) 

CwLL  SYMBOL  (  .5 ,9  .b  ,.l  ,  HEAD,  0,7  01 
C.AU  1 YH9C  L  (  .  6 , 9  .  E  *  .  1  *  HE  AO  i  ,  3  ,  7o  > 

CALL  f  Y  M  9  O  L  (  •  '•  ,9.9  , .  1 ,  Ht  AO  2,0,  70) 

CALL  AX  IE  (  v.,  ,L.  , 

c  5HGK  CUP  CElAY/CM  (♦>  (NSEC), 25, >X  *  Y  C  .  » YM  IN  ,  0Y  *  Y  Y  > 


CALL 


CALL  AX1M  C.  ,  j,  ,  lSH^fcOUENCY  (MHZ), -IF,  2  *t.« 

CALL  LINE<  FM,VGM,J1,1, 0,1,XMIN  ,0Y ,YM1N,CY  ,  ,0S ) 
CALL  FLn7(i-j,,Ll»»-2) 

CALL  SYMPOLf.P  ,9.0  , .  1  ,  HEAD  *  j  ,  7  u) 

CALL  SYM9UL  (  .9,9.6  ,.l  ,  HE  AD  1 , 0 , 7-J  > 
f.-L  L  CYM9QL  ( .5 ,9.4,.l  ,  HE  AO  2 ,3,  ?u> 

CALL  AXIM  0.  ,J.» 

I  2LHG9CUP  DFIAY/CM  (- )  (^^LC),^5,X 


,  L.  ,XMI  N  ,CX  ,YY) 


CALL  AXIS ( L.  , 
CALL  LINf  (  FP, 
YMIN=C . 

0Y=  cd. 

CC  27  1=1,  II 
IF  (Ff(I)  .  >  T 
j. r(F.p(j>  . r.r. 


lEiM-Pfc  QLENCY  (  MHZ  >  ♦ 


(M  JLC  ) 
-1‘, 2 


,11,1,  0,1,XMIN,UX,VMH,C  Y,  .ue  ) 


30  0.)  KP(I)=iliL. 
2  l  Co.)  hP(I)=20lw. 


25, XX  , EC.* YMiN.DY.YY) 
,0. *XMIN,CX,YY) 


CO  l  liw- 

CML  FLOY ( lo., G.,-3> 

L  *- 1  l  AXIS  G.  ,27HRAO.  RES.,PINUS  WVE  (CPP  S )  ♦  2  7 , 10  .  ,  90  .  , 

CuVui-u:;,.  .lM-F^tQUtNC  Y  <PH7>  ,-lE,7  . 1 .  ,XMI  N  ,  CX  ,  YY) 

I  at  l  J  N!  (  IP  ,  <i'  ,  I  1  ,1 ,  u  ,  1  ,X  M  3  N  ,  LY  ,  YHIK  ,L  Y  ,  .  jt> ) 


!.  L 

?*5  J-l. 

Jt 

’  f 

(  'P ( J 1 

.r,T. 

?  0  j  . ) 

AM(  Jl  =3  j  , 

• 

IF 

(RP ( J  ) 

.  L  T  • 

2  C  0  C  .  I 

RP ( J ) =2  ! 

00. 

IF 

l  ’ll  ( J  ) 

»•» T  • 

SCO.) 

R  T  ( J )  =  3  0  0 

• 

lw  • 

1  F 

< r  T  ( J  > 

.GO. 

seen.) 

RTU>=2. 

LCff  1  Nl)E 

l-'LL  FLC1  (  1/  .,  b  .  ,-3> 

e-LL  «1L  (J.,  0.,?7HRan.  RES.,  PLUS  PAVE  ( CHF  SI  ♦  2  7 , 1 L  .  ,  90  .  , 
YMJN.OY. 1j.» 

CALL  A XI f  (  C.  ,0  •  *  1  E  FFR  t  QUEN  CY  (PHZ),-1S,2  , C. • XMI I ,CX , YY> 

fill  l  L  INC  <  F*,RM,  Jl  ,1,Q  ,1,XMIN,  CX,  YHIN.OY  ,.031 
i  ALL  F  LU  T ( 17. , D.  ,  -  3  > 

CHL  AXIS  <  0.  ,0  .  ,221-RAU  .  RES.  TOTAL  ( C  FMS  )  , 22*10. ,SC.,YMIN,0Y, 10  .) 
CALL  AXIS' j  .  »lpHFRfc  QLENCY  (PH7),-1!,2  ,0.  ,XMIK  ,0X  ,YY) 

(III  LINL(F0,-;T,I1,1,  (J,l,XFIN,CX,YMIKfuY,.i)b) 

1  1=  p 
j;  =  p 

YPlNs-lwlJL. 

YMN  =  -25w  . 


•;  t  =  _  j .  j . 

c  y  = ; . , 

.C  "2  1-1,11 

lr  IFX(I)  ,lT.  - 1  G  C  UJ  •  I  PX  (  I)  =  -10  J00  . 

Ir  (PX(I)  ,LT.  -2S3.)  PX(I)=-25u. 

U  (Pxil)  ,LT.  1CCCG.)  PX  C I )  =  1  c  3  J  J  . 
ir  irxn  . >■; t .  2E3.i  Pxti)=2r,G. 

Ct.N  I  INUl 

rxi  plot  <  l?.,  l  . ,  - 

v.  ►  L  L  AXIS  (u.,3.,22  FPA  U  .  REAL  TOTAL  t  C  FM  S  )  ,22*10. ,*:(•, YMIN«DY  ,10.  I 
uUl  A  Y  If.  (  c  »  ,  j  •  ,  1 5PF  Rt  QUENCY  (PHZ),-15,Z  ,0.  ,X MI N  ,  CX , YY) 

Cflli  l INI  (  F°,»X,I1  ,1, C  ,1,XM1N, OX, YH1N,D Y ,. L8) 
rC  -.s  1  =  1.11 

IF  CE^P(I)  ,lT.  -o'!.  )  SERP  (  II  =  -bu  . 

IF  <SSRP(I)  .LI.  -H).  I  SFkM  ( II  =  •8L‘  • 
l*i.NT  INU:. 

YMK  =  -t:  . 

r  y=iu. 

C  iLL  FLIJ1  <  1*.  ,  L  .  ,-3) 
f  -LL  :  Y  M  3  O  L  (  .  ‘j  ,  v  .  f*  ,.l,HtAO,u,7C) 
f  -»LL  SYMBOL  (  .  i  ,9  .c  ,.l  ,  HEAD1 ,  J,  70  I 
fALc  S  Y  M  *>0  L  (  •  ‘  ,9.  A  ,.l  ,  HEA02  ,J  »  ?l  ) 

C  ‘'LL  AXIMC.,  .,261— INS.  L 0 b S  ,  PI NUS  fcOVE  (C3),2b,  t.,S0.,YMIN, 

<  •'  Y  ,  it. .  ) 

(  'LL  JX15  <  c.  ,.i  ..lEPUcOUtNCY  (PHZ),-15,2  ,  0  .  ,XMI  h  ,CX  ,  YY  I 

i  '  ll  LlNf  (  FP,SCFP,Ii,l,U,l,XHlN,OX,YPlN,CY,.Ot) 

CALL  HOT(  in.,  0.  ,-3> 

tall  SYMBOL  (  ,  i>  *  f . o  ,.l  ,  HEAD  ,0,  7  t) 

(ill  TYMcOL  (  .  >,9.fc  ,.l  ,  HfcAOl  ,J,  ?{  I 
i'll  S  V  M 17  O  l  (  •  j  ,y  .A  ,.l,  FilAu2  ,  J,  70) 

L7L  i  AXIL  <  L.,0  .,260-1*5.  LCbS,  °l  US  HAVE  lCc),2t>,  t  .  ,<s0  . ,  YMIN, 

!  Y  ,  *  L  . ) 

I-U  »-Xlb(  t .  ,< .  ,  1SMF  RtQLENCY  «PHZ),-1E,2  ♦  C .  ,XMI  K  ,CX  ,YY) 

Mil  l  INE(  FM.bEPP,  11, 1  »  0,i,XMIN,nX,YMN  ,CY,  .CF) 

(Alt  F  MDPL  T 
L  IOP 


F  I  M  7  '• 

l  •  o  r- 

FC«?PA  1  { 1  HI  ) 

F  C  ~  f  I  l  ( :>  X  ,  M  H- 


F  l  P  3  I  C//ii  (b 
r  L!;PAT  (///F  IX, 
FT- PAT (///  b  iX, 
r ( f  ri  \  (///  (1JX 
FLFKATC/// <i 
h  r-vF'tr  (  //  duY 

F  f  a  pa:  (///  dux 

flrpmt  r*i  i 

FCPMAT(//(  1  ,  . 
F  C  ->  P  A  T  (//  (  lax, 
FCRPA1  C*  «-s 


■,ri9.7/5X,"  T  1  =  ”  ,  E 1  £  .  7  / 
5  X  L  =  ** ,  E 1 5  •  7/ 


,E1£ .7/ 

E>  ,**POOE  = 


■,F3.0) 


Pi 
i  2 

f  7 

£  5 
1  0  * 


F  F  :  MAT  (//  (  1  jX,  ”1  1  (*• 
FC-  t'Ml  <//(  1-iX,  "  AC 
r  C  F  f'  AT  (✓/(  i.:X,"  PC 
FC‘ f'AT  (/  (IlX,"  F  =  *• 
F ( ‘  MAT  < 7  A 1 b  » 

?  N  r: 


,l4d  =  M.C15.7/»» 

1 1  ♦••)=  "fF  15.  ?/ )  ) 

e£lb. 7, 1JX, "GROUP  tFL/Y  (4) 


*  ,E15.7/  >  ) 


FUNCTION  COTH(CA) 
CCTF  =  l  ,/TaNH  (GA) 

F  I T  L  t  N 
F.  N  f; 


FUNCTION  CCT  (C  A  I 
C  C  T  =  3  .  /  T  AN  {  C  A 1 
=  £1 UFA 

fn: 


FUNCTION  FllCft) 

COMMON  D»EC«T1«j»C»UYT  »PI*ENM.FCA 

A=2.-CA*7l 

CAISCA) 

IF (A  ,G1  .  1 75. I  GO  10  5 

Ai=  CCIHC  CA*tl)+l  .-eXPC-c.*CA*Tl»*  (CCTMCA*  £  l ) -1.  > 

GC  1C  0 

£  A  1=  c  C  r  H  (  ( 

e  CONTINUE 

3 1  =  - 1 ,  4C0TH(j:A*£L  c  «C*UYY)**2 
IF  (A  .m.  6 /  i  •  )  CO  10  2 

C  1  =  t  X  P  (-2.  *C  A*  t  1)*  (-1. 
v  -CCTMl  A  4  c  l  >  *  (C*UYY  )♦♦?) 

GC  10  3 
2  0  1=1. 

3  fCNUNUf 

Fi=  ATAN(C*UYY*At/  (61-Clll  /C/C 

Ir  (FT  .Li.  j.)  F  7  a  F  T  +  A  3S  ( P  1/C  /  0  > 

FT=FT  +  A3G  (ENM*  PI  /U/U  ) 

F  C  A  =  i1  *C*U  YY*t'0  1  (C«CA*u)  -31+Ci 
F£TU=N 

fn 


Fi.Nv.lIGN  l3v.H(CA) 

CGCFaC. 

IF  I C  A  .lc.  ?’’■*  L  •  1  CSCH=1./SINH  CCA  > 

*  T  (. c  N 

F.N3 


FUNCTION  CINC(CA) 

Fl=  1.14159265- 
C  ]  M  C  =  ( C I N  (  F 1  #C  A  )  )/  CPI’CA) 
F  T  L;  N 
'NO 
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Reproduced  from 
best  eveilebte  copy. 


KiKCIIGN  f--7j.(Cfi> 

CCHt'CN  Jt  E  I,  11  ,2  »C  ,Ut  1  ,  PI,  cNM,  FCA  ,OPN 

xnr'|((U;‘ai':  )’n-C*n^LMfS'.HCA*£l  ))  **2>*S1^I(C«CA*D>^• 

X  (  L*  ( C  *  *  2  ♦  i  . »  *UOl  F  <CA+  1 1 )  *c.L-  (uSCH  (CA*  El  )  >*♦£  >  *C*C  CS  (C*CA*i)>  ) 
x '  ■.  xp  <-ca *  )) 

</  ^  Hi.  )  ♦  C*'Ol  H  (CA*  t.L  >  ♦  EL  *  (CSf  F  ( 1 1-  *1:1  *  )**2  »*iIN  i  U*i  A»0>  ♦ 

<  <,  •  <r  +  *2«-l  .)  ♦  El  *  <ISCH(  CA*EL  )  )  *  *?)  *C*CCi>  <C*CA*C»  > 

X*„-  t  XPI-PA*0) 

?-Ac.Z  v:.t <J  +  ?  .*71  )) 

li-  <«!  .  GT .  o7  .  )  GC  TC  S 
IE  (GFN  .  r.  G .  .  .  > 

<F1".  =  (L*~2*C0TH  (CA*EL  >♦!.•  I  -  el  *(CSCH(CA  *El)  )  «  *2  -  (  Z.Ml+O)  ♦ 

X  <CC!H<(  A*_L)-1.) 

Tf  <C*N  .^C.  1.) 

x  F13=-(D+2.*ri)*-  (CCTHIC  A*£L>  -1  .  >-EL*  (C  St  F  <C  A*E  l  >  )  *  *  l  - 
X  0MC**2+COTH  <GA«  r.L)> 

XFT-  =  (2*?J*T1>:,J<C**2*GUTH(CA*El>+1»>+C**2',£1*(  CSCH  «CA*EL))**2  + 

x  r '**z*a+  <  corn  (Cam.) -l.  > 

If  <  C  F  N  ,lC.  i.) 

y  F  !  -  =  -  ( U ♦ 2 • * T 1 ) *  ( C*  *2+U OTH ( C A*c  L> ) -EL  *  1C 5 C M C A *£ L )  )**2  + 

r  i  i  =  f  ■»  -  \  P(  -CA‘-  (C*2  .*  T  1  )  )#(C,C02(C'*CA-fC)  *F72-SIN(C4IA*C)*FT‘+)  /G 
x  +  FT  C / C 
fC  TC  J 
F  .  1  =  V 7  2 1 C 
<.  CN  i  INUt 
•  ».  ri'N 

i  Nr- 

f  L  N  C  T  i  UN  Gf>  (i,A) 

CCMFLFX  f.ii 

,  ccwvrx*  n,ei  ,T.  fC  'UYY  •PltEKPt  FCA  .OP*  fE  If  (EfcvELl  (3D  tPE(3d)  ,  AA  t  30) 

N  =  •  M 

r  if.  ^LX  («;  .  .1 

i.  £  -  CMP  L  X  ( C  OS  <CA4  I*PE<  I  >>  ,-SlN(CA*I*PE  <I>  >) 

«  -r  *siNi: (. 5- a..  <:>*ca/pi)*eta*’i*sqrt<eii mi*cs 

L  CNT IXUF 
GfYaCA3£(-» 

Pi  TUP N 


f  lNCTICN  Gc  (Crt  ) 

rcvrCN'  0  « E  L  *  T 1  ,S»C.»'JYY  ,PI,£NK,FCA 
i"  -  sAri  (GAY  (CM  ♦cXP(-  U  A  ‘  L  l  /  F  Tl  (  CA  )  ) 
S  1  LFN 
f  Nil 
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fL^-'OL'iJN.  MT!#-N'  (f>  «  X  «  N  «  F  '  l.  G  «  F  r  Nli ) 
l  i  (*  i.  M:  I  C  N  -  (3  )  ,  y  (  3  ) 

PI= 

FL-i-lf  fj-ra  G  W ( N-l  ) 

F  =  F  3EG  ♦  . E  v  F  3  £'. 

IM’:Fnrj  It),  c) 

M  =  mim:  -i 
FK1=N-1 
MMI  =  NI-? 
re  :i  i=l,K'll 
*(!)=., 

Tr-  (T  •  „  0  .  1>  Fy=l’.*'«<l)^b.*;|2)-R(i)|/t. 

if-  (  T  't->,  .'IMD  X  =  (  -  •*  <  N-2  )  ♦  fc  •  *R  (  NH J.  )  43  ,  «  £  (NH/6. 

if  (1  .-If.,  1  .  r f- .  I  .to.  M M 1 )  GO  TO  2  3 

<-«  ccmimk 
F  1  =  F  '  0 

r.i.  J  IJ-i,NI'U,2 

IIF4tl-iXI/((  F  1+F  )£L»  **  <-F**2 ) 

X  f  2 .  '  O.  (IF  )  -  k  X  )  /  (  rl  «*2-F**2) 

F i=r 1*2. *F C-L 
CCMTNtt 
f  _  N  -  F  E  N  ij 

1M1NC  .tO.  U  )  r  r  N=F  t  NC-FC  EL 
Mi  1  =X  ( I )  ♦  <  ( M 1  >-•=>)  ✓  «F£N**i  -F  **2  » 

X  -  (  ■*  ( *  )*«»/  (  FHo**  2-F"*2) 

>  tl)  =  FCa/i.  +  <  (I) 

if  (INC  .  tv),,.)  GO  TO  ij 

>  (  T  )  =  /  ( I )  4  .  ^  OF  L*  « ( ^  (  NI  >-(?>>  /  (FF  M**  F-F**  2) 

X  ♦  U(N)-FX  )/<FFNC*»2-F«*2  »  > 

o  .  X  'T)=2./PI-F‘X  (1)-»CX/H1*AL0G 

X  Ul.-F/KNi)/  (i.tF/^tNO*  (F  +  F5CG  )  /  (F  -F  EcC)  ) 
MF+FiEt 

33  CCNTINUi 

X  (1C.  -■  >  (  1)  -A  {  .  •»  X  (2)  ♦  3.*X  (  3)  >  /o. 
X2=(3.4XllT^-o.4X(2l-X  (3) )  /<,. 

EL  21  Is  2  fNH> 

X  Ts  (-X  (1-2  )>  3.  ♦>  ( I-i»  +0.*X(  i»  -XII 4-1  )  >  /*fc  . 

X  (]  -2 ) - X 1 
Y1=X2 
XT 

3*  CCNTTNUC 

xiN)s(l5«->(N>1i)-lC.*X(hf»2)  +  3.*X(N-3n/e. 
x  (N-l)  s  (  ?.  INM1ME/X  (M2  l-X  (  N-x  )  »  /  c  . 

/  (N-2)  =  X2 
X  ( M  -  3  )  =  X  i 
;  T  L  -  N 
‘  N  t; 


25 


B.  Dispersion  Relation  -  Generalized  Volume  Waves 
A  computer  program  has  been  implemented  on  the  CDC 
6600  at  Hanscom  AFB,  MA  which  obtains  the  dispersion  rela¬ 
tion  and  group  delay  for  volume  waves  of  arbitrary  orienta¬ 
tion  of  the  biasing  field.  Plots  on  the  Calcomp  plotter  at 
Hanscom  AFB  as  well  as  computer  printouts  are  produced  for 

the  wave  number  and  group  delay  as  functions  of  frequency. 
One  can  choose  the  solution  mode  for  the  dispersion  rela¬ 
tion  (42) .  The  frequency  range  is  determined  by  the  pro¬ 
gram  from  (43) . 

The  use  of  the  input  cards  to  this  program  now 
follows 

Card  1-  Hq,  t^ ,  d,  0,  l,  <J>,  mode  number 

Seven  quantities,  separated  by  commas,  are  supplied 
here.  All  lengths  are  measured  in  meters.  Angles  are 
measured  in  degrees.  Their  orientations  can  be  determined 
from  figure  1.  For  the  fundamental  mode  the  mode  number 
is  o.  Columns  1-72  may  be  used 
Card  2-  heading  for  plots 
Card  3-  heading  for  plots 

Here,  there  are  two  lines  supplied,  in  order,  which  are 


used  as  headings  on  several  plot  frames.  Columns  1-70  may 
be  used  on  each  card. 


C  Th--;a=-..  uE.0  -  P1CKWAR0S  VClLPE  FAVEi 

C  PH  1=90  OcG  . . .  FOPVA"'C  VCLl ►£  WAVcS 

PROGRAM  VO  L'1(i'  NFIT  *011  IPUT) 

C1M2KEIGN  F-*CGIC  (2  ) 

CIPrMSIGN  F  (  i.  i  1  >  .FP(  12  1)  .FM  (12  )  ,  C  fP  (  1 1 1 1 )  ,CA  M  <3  t  0  1 ) 

X  V  C  F  ( 1 2  l).WHH:,i)fElAP(;2  l).rLAP<.  2  1) 

C1M.NFI0N  F  A ; 3 < 7 ) ,H5AD2<7 > 

c;  C  ^  h  CN  0,a,  T-  ,S  ,C  ,LYY  ,L)YXYYI  ,  FI,  2NP,c(.  A 

F  Ar  *  ,H,  T:  ,  Q,  7H£  I  ,t  L  ,  PHI,  PPM 

;3:AP  1  ;j  ,  H  r  '  0' 

'  M  1  C  C  *  H  f  .  j  •’ 

F:IM  b  j 

P  >  7  NT  6i«H  ,  .  j.,  (.•  THEY  ,  t^PHl  «£Kf 
F'-I  M  f, 

f.  FT  F-OGIL/ 'HSi  lPA-<£S  M  J  .  W  -  IN  £  ;  r  G  / 

IHl  FLT1  j  -  IP-  Of, TO  ,U  J  .  ,12.  ,1  ) 

P I  =  2 . 1 A 1 1 .  c  P 
THLT  =  th  ",  FT /i  <  (? . 
r-HsPH-'PT/i*  . 

T  «F  r “ L 
c  l.=  1 

V  j.r  MF 

L?-  . 

GAM -  2. F 

FLO-'.c  SOP.  (H*  (hH?V  .  (SINChT))  2  CIMFH))  i  >  > 
,0'-  sop;  (**  ( F  +  1 7b  (.  , )  ) 

F -  «  =2. 

f  '~r-  =  A  1*4 ~  <  Fl.C)  +25. 

NF=  I  NT  C  P  HI  )-I,n  IFLf.l-ij 

M=NF/2 

in  i=i,  m 

U  ‘  F  (1  )=FO(-.G+  CI-l  )•  FUEL 

IF  <  F  ( 1 )  .Li.  FLC)  F°  1  N1  *  ,  **F*  tOU  ENC  1 ICC  LCW" 

IF  <  F  ( NF  )  .GT.  F  H  i  )  P«IN7  +  ,  T  FLO  UENC  !r  E  7CC  UGH" 
U=GAv**2*p  £>1 
K  =  i 


f:  F  =  r  (K) 

T  f=  <GAM-‘H)  1  ?-EF  **?. 

T  c- GAP+EH* i - 

UV>  =TU/TH-  ( ‘  IM  (TFcD)  - ■  * 2  +  1. 

LXX=TL/TD*  USj  MCTH5I)  “SIN  (PHI)  >  '  2  + (CCS  <  IF 
U*YYX  =  -2.«-  IU/,'0-  CCSUFETHCOS  (PHI  J-SIM  UE 
LUX  YYI  =  T2/7  <•'  J  I  N  ( 1  HLI  ) ' SIN ( PH  T ) 


T)  H*2>41. 
) 


/('  YY 


T  *-H  =  UXYYX*  #  -  .  M  X  X®  U  Y  Y 

L  =  2 

IF  UFM  ,LL.  •  )  F c  1 N  I  - 
IF  OLH  ,L  t .  .  )  GC  I  0  33 
C  =SCPT  (TEM  )  ,  /(i YY 
L  =1 

IF  (L  .!)  GC  1C  ? 

L-l  . 

GC  TP  3 
?  =- 1  • 

G  C  M  ">  I  MJ  ?: 

IF  (J  ,f  I  .  1!  C-C  TU  5  3 
CCO  =  ATAMC,  C  UY  Y/  (U;*UYY)  *' 

ip  (C a o  .lt,  .  )  caomao+ab: 
OOMAO+A3E  (  5  IK*  PI  n./M  ) 
CCNTTKU- 

xF  (I  . wP.  ' )  GO  10  11 
IF  (J  .lT.  *)  GO  10  51 
1 F  (  L  .EU.  1)  AO=CAP(  1-1) 
IF(L  .EQ.  c)  0  A  0  =  C  A  M  (  J-i  ) 

P  =  1 

l  f  K  ■  '  I  Mi 
CC',  =FT  (GAO  ) 

Pi  I  M  •'.tAr.C-.I.FLA 
IF  IAES<  (C  Fi-'AC>/CAO>  .LT. 
C  A 0  =  C  A  1 
H  =  H  +  1 

I  P ( P  .GT.1F)  GO  TO  35 

r  r.  i  o  5 

IF  (L  .  rO.  O  GO  7C  2  ! 

L*i  =  CAi 

IF  <A3S(FCA  )  .Cl.  1.1  GC  Tl 


N  a  C-A  T  I V  E  SOI  ARE  R(OT‘ 


!  C-C  TU  5  3 

UY  Y/  t(<;*UYY>  *-  2  +  LYXYYI  ' 

.  )  c:au=lao+A3s(pt/c/d) 

IK*  PI  A,/U  ) 

GO  10  U 
x)  GO  10  51 
0 AO=C  API  1-1) 

•; ao=c am(  j-i  > 


J  ,FC  A 


l.))/C/0 
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j.  )  GO  TO  35 


IF  (CP  .W. 

FFUIpFF 
C  pp  (I  )=CA 
i -r  ♦! 

FklM  '  ♦  "C  A-  **  «  IF 

L  =  2 

rn  tp  o! 

C  «  =  C  A 1 

IF  (AFStFCA  )  .FT.  1.1  GC  TO  .<5 
IF  (CA  .LT .  0. )  GO  10  35 
F  ¥  (  J>=EE 
i  AM  (J)sCA 
J  =  J  +1 

FAINT  *,mCA=  ".CA.-Fr  " ,  EF 
IF  ( <J  .FQ.  I)  GC  10  15 
If  (J  .GT.  I)  GO  7  C  31 
I=J-1 
K  -  K  - 1. 

J-J-l 

K  ■  K  4  i 

IF  (K  ,LF.  NF)  GC  10  50 
F  < I N T  60 
1 1- 1-1 

GC  T0  Z+ 

P’.JM  *."I1 .3-1 ICft  00  s  S  NOT  CONVG  \GE.F=  ** ,  E  F  »  *  S  =  *  , 
7  f  (L  .  tO,  J  1  GO  TC  15 
L  =  ? 

GC  '0  2 
i.  CN1  IN  Ur 

D  ■  1  NT  6j,  (F=MlL  >  *CAP  !Ul,LL=liL  ,  .> ) 

F'lM  64  ,  (  FM(Ll  )  *CAM(  LL)  ,LL  =  1,  Jl,  •’) 

F  ■  1  M  o  j 
lC  11  1=1,11 

K0C  (I>=2  .  P  /CAP  (1)  *  l.Eo 
GC  12  J=l, J1 

IE  IAMJ)  =  ?.*”!/(  AM  jri  ,tu 

MM  65  ,  (  FP  (i.  L>  ,  E  LAP  I LLI , Ll  n  ,11 ,5) 

FFIM  fcb,  (FM(Ll  )  LAM  (  IL )  *  L  l  =1  »J1  »  5 1 

C  C  <-1  1  =  1,1. 

IF  (I  .L3.  .)  VGP(I)-5./PI*(CAp(2>-CAF(  m/FCEL 
IF  (I  .EG.  IU  VGF  (1)  =5./PI  (cAp(  ID-CAF  (Il-l  ))/FCtL 
IF  (1  .  NF_.  _  .  ANC  ♦  1  .Nt.  ID  VGP(I)= 

<’>  */Pl*  (C*»P  (i*  )  -CA  F(  I- 1) )/  FOEl  ‘.5 
CONTINUE 

F'lM  ft  7  .<FP(1  >,  VGF  (I  )  ,1*1,11,51 
CC  c2  J= 1 , J1 

IF  (J  .EG.  D  VGMJ)  =  5./PI*  (CAN(2  )-CAN(  1  H/FC  el 
IF  (j  .dO.  J 1 )  VGM(J)s5./pi-  (iAM(  Jll-CAIMJl-l  n/FCEL 

IF  (J  , NF ,  :  .  A  N  C  •  J  .Nt.  Ji)  V  GM ( J )  = 

(  /PI*  (CA M  ( J4 1  )  *C  AM  ( J  -1 )  >  / F C .  L '  .  .> 

CONTINUE 

F  ;  I  NT  o5.(FM(Jl  »VC-K(J) 

VN1N=2300 , 

XNTN=s2  . 
xmn=7j:<:.  . 

LX*  i  C  t  • 

c^r  ?cn , 

V  I”  T  N  =  . 

r.  i  = r.  ••  i  j.'i 


I  N  T 
FFINT 
CG  <-1. 
IF  (I 
IF  (I 
IF  (I 


:  v*  i  .-'.6 

i:c  13  i=i,  it 

IF  (CAP  (I )  ,G  ,  1.E6)  CAP(  I )  =  *  ifo 
IF  (FLAP  (I  )  .Gl.  1(10.  )  FtAMi)=10CC 

VCF(I>  =  A3S  (VG‘MI  )) 

IF  (VGF(I)  .  G*  •  l.lj.l  VGPU>  =  :  , 

rc  s  j=i,j: 

IF  <C;M(J)  .C  .  1.E6)  CAM(J)=.,-6 

IF  (lLAMU)  .1,1.  i  j  .  ?  .  )  El  AM  (J)  =1  Ci  „  . 

VGM  (J>=4?$  <  7  GM  (  J ) ) 

If  (VGM(J)  .  G.  i;M.)  V  GN  ( J  )  =  1  jj., 
7=11. 

X  -  i  . 

>'•  =  !!;. 

X  X  =  1 G , 

CHI  FLOTC  1.5,  0  .,-3» 

CALL  SVf»30L(.-  #  r  .6  ,.l  ,  HEAQi  ,0,7  > 


C  £  L  L  SYMROK.  , 9 .4 , .1 , He A02 , j , 7  ) 

C  A  L  L  A  XIL  (  l  •  *  .,2cFWA  J  E  NUMBl-  K<4)  <  1  /  M  f.  1 ER  >  ,2  o  ,  X  j  .  ,  YM  IN,  DY  ,  Y  Y  ) 

CALL  SXii(  l.  ,  .  , laFFBBQUSNCY  (  MH7  )  ,  -  1 :  ,  2  ,  C  .  ,XMI  b  ,cX  ,  YY) 

C «L  L  LINf  (fM.CAM, 

X  Ji  ,1,  U  ,  1,XMIN,  CX,  YMTN  ,,,Y  ,  .  Qc) 

f ALL  FIC1  (,o.  , „  . 

CALL  SYMBOL  (..  ,S.fc  ,.  1  ,HcADl  , ,  7  .  ) 

CALL  SYMBOLU  , ;  .4  ,.l,HEAD2,  J, 7  ) 

CALL  AXIE(  u.  ,  .  .  ,2bFWAVE  NUMBER  K(->  <  1/  M  E  7  ER )  .26 ,  X  0  .  ,  YM  IN  ,  DY  .  Y  Y  ) 

CALL  AXISU.,  •  » 15 HFRfc QUENCY  (MHZ),-1£,2  tC  .  ,XMI  F  ,0X  ,  Y  Y) 

CALL  LINf  (  FR  »>;  AP  « 

X  Il.i,U,l.XHiN*C**YNIN,C.Y,.04> 

C  Y  =  *  ,• . 

CALL  FLCT  ('3  ,  L.,-4> 

CALL  5YM30  L  f  »  -  »c  »b  , .  1  ,  HEAD1 ,  J  »  ?  i ) 

CALL  SYMBOL  (.  j  ,9 .4  ,.l ,  HEA02,  rt  ,  7.  ) 

CALL  AXIS  (  ...  ,  ,  LAHUAVc  LENGTH  (4)  (MIC^CNS)  ,  2A  ,XX,Ew»fYMIN»DY,YY  ) 

CALL  AXIS  (  *  .  ,  .  ,  15HF  RtQUENCY  (  MHZ  )  ,  -  1  ;  ,  2  ,  C  .  ,X  MI  N  ,  CX  ,  YY  ) 

CALL  LINE  (FM,t  LAM  ,  Ji  ,1,  j,  1  ,XMIN,  OX  ,YM  N  ,CY,.  Cd> 

C-'LL  PL  CT  ( L  3  ,  0.,-3t 

CALI  JXISI,.,  , 2  hH  XA  Vc  LtNG  TH  ( -)  <  M  ]C  5  CN  S  )  ,  2k  .X  X ,  E  t . ,  YM  IN  ,  3  Y  ,  Y  Y  > 

CALL  A  X  IS  (  t  •  »  . ,ibHFRtQUENCY  (MHZ). -15, 2  ,  0  .  ,XMI  N  , CX  ,  YY  ) 

CALL  LINE.  (FP,  EL  AF  ,11 , 1,J,  1.XM1N,  JX.VMJN.lY,.  C6) 

CALL  FLOT(  It  ,,  0.  ,-2) 

CALL  SYMPOL  <  .  ,  3 . 6  1  ,  HE  ADI ,  J  ,  ?u  ) 

CALI  SYMBOL  (  .>  ,  9.  *.l  ,HtAU2  ....  ,  7.  ) 

C*LL  A  XI S  (  ^  , 

X  H  G  «  CU  P  DELAY/CM  (♦)  (NS  C  )  » 2  5  ,  >  X  »  b  C • »  YM1N , DY , YY ) 

CALL  A  XI S  (  ‘ .  ,  .  , 1 BHF  R  fc  QUENCY  (MHZ), -15, Z  ,  ti  .  ,XM  I  M  ,  [  X  ,  Y  V ) 

CALI  LINC(FM,/GM,Jl,l,0,l,XMIN,DX,YMn,cY,.Ce  ) 

CALL  FLCT  C  1  »  C  . ,  -  3) 

C«LL  SYMBOL  (  .v  ,b  ,.l,  HEAQi  ,  i »  72  ) 

CALI  SYMBOL!.  *  <•  .4  1 1  HLA02  ,  J  ,  73) 

CALL  A  X  T  S ( C  .  ,  .  , 

X  L’sFGRGUP  DELAY/CH  {-)  (NIFC),25,)X  ,  %  L  . ,  YM  IN  ,  0  Y  ,  Y  Y  > 

CALL  AX  IS  (  i  ,  •  ,  1 b  HFRt  QUENCY  (MHZ),-1-,Z  •  (l .  ,X  MI  N  ,  CX  ,  /  Y) 

CALL  L  INC  (FP,VGP  ,11,1  ♦  0  ,1  ♦  XMIN  ,OX  ,  YM  I K  C  Y  ,  .  l)  6  ) 

CALL  CNDPL1 
S  1  OF 

fc  F  C  7  MAI  ( 1H1  ) 

€1  F  C  F  MAT  (S  X  ,  H-  ”  ,  t  lb  .  /  /bX  ,  ”  T  =  "  ,  E 1  E  .  /;>,**  L~  ",El5.7/ 

X"  TFETA=",F.5  ?/5X,"  L=  ” »  ei 3 .  2/5  X  ,  ”F  F"  =  ” ,E1 £• 7/ EX  ,”MOD£-  ”,F6.0) 
fc2  F  CRMA  T  (/  /  1 1  <  5  - 1  E  , 5  /)  ) 

bo  FCFMA1  (///£  :X,*"S=1*V/  (  10X,**F  -  **,E15»7  ,1C>,"K  (-)=  *,,tl‘3.Z/)> 

FC  FMAK///5  X.  “Ss-r*//  (1  3X,MF=  ”,  £15  .  / »  1 1.  X  ,**K  (♦)  =  "tti.a.7/)) 

ti.  Ft  7  MAI  (///£  X,”S  si-//  ( 1.  X,”F=  ",r.L5.  1  ,loX,"LAM-)r  ",E15.7/n 
rl;  FC^MAT  (///=  X.”G  =  -1”//  (13X,”F=  ”,  E13  .  ? ,  1  i.  X  ,”L  AM  (  4  )  r  ”,E13.7/)) 

S  7  F  C 0  M  A 1  ( ”  f=  ”  ,  F  lb .  /  ,  ”  G  <  CUP  Dr  l  AY  (-)  =  ”,(1E.7) 

F3  f  C  P  M  A  V  ( /  ( 1 ..  X  ,  ’*  F=  *,fcl3,7,l-X, "GROUP  (t  LAY  (4)  =  *',Elb.?/>> 
lo.  FORMAT  (7A1  c) 
f  t\n 
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nun  Mi* 


FLNCTION  COCH(CA) 

COMMON  D.FL.T,  ,  C  ,C  .UY  Y  .UYXYYI 
IF  ICA  .GT  .  1.3)  GO  10  3 
IF  <CA  .IT  .  -17.3)  SO  TC  s> 
CCTH  =1  ./ TINH  (CA) 

GC  7 C  2 
CCTFsl. 

GC  TO  2 
CrTH=-l. 

CONTINUE 
x.T 
F  NO 


FlNlTICN  COM  A  ) 

CCMF'ON  JjiL.Ti  ,S,C*Ur»  ,LVXYYI 
CCT  =  1 .  / 1  AN  ( C  £  ) 

Fr  TC^N 
t  NO 


F  C  NC  T I CN  FTlCt  1 

r. r.MKCN  O.'IUT,  ,G.C .UY Y ,UYXYYI,PT, SNW.FCA 
A  *  ?  v  C  A  7  l 
A-AGS (m) 

I  F  1  A  .GT.  C7s.  )  GO  TO  5 

ai=  CCTHU  A*EL»*l.-EXPl-?.' CA-711  ICCT  H  (C  A*  F  l ) -l ,  ) 


I  F  1  A  .GT.  £,/•?.  )  GO  TO  5 

rd=iC6  CPTHU  A»El»M.-EXPl-?.J  CA*T1I  ICCTFMC 

M=  CO.  H(f  A*r  U  +  i. 

C  CM!  IMJF 

M*  (CCi  H<  C  A*il)-l,)*S4  UYaYYI-1 

fl=el  ♦CQThCCA'  si  \*  (  I  C+UYY)  *  2+UYXYY]  * 2 ) 

JF  «A  .Gl.  6<5.)  GO  102 

C  1  =  r  XP  (-2.  '  CA  TDMIfJOTHICAIUH.)^  UVxYYI-1. 
<  -'OTMCA-  .U'  t<C*UY»)**2  +  0Y>YYI  “ZM 
G  C  T  0  3 
Cl=  . 

(  cm  ikuf: 

F]=  ATANILMJYY*  Al/  H1-C1M  /C/[ 

ICT<U  S W J : r-  a>  P  =pt  ♦A3L(Pl/C/0> 

FT  =  F;  ♦  A  3  S  <  EMM'  PI/C/C  ) 

FCA  =  Al-»C*UYV„OT  1C  ’CA*  01-31+C  . 

F  t  1  L^F1 
r  NO 
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C.  Hyperbolic  Current  Distribution  -  Surface  Waves 
There  is  a  surface  wave  program  operational  on  the  CDC 
6600  at  Hanscom  AFB  which  employs  the  hyperbolic  current 
distribution  described  earlier.  The  program  produces  plots 
and  printout  for  quantities  of  interest  for  surface  waves. 

The  option  for  a  truncated  array  of  normal  modes  is  eliminated 
from  this  program.  Otherwise,  the  program  is  similar  to 
the  basic  theory  surface  wave  program  described  in  the  Interim 
Technical  Report.  The  number  of  input  cards  is  the  same  as 
in  the  earlier  program  and  the  use  of  most  of  them  is  the 
same.  We  here  describe  only  those  cards  whose  functions 
differ  from  the  earlier  version. 

Card  5  -  AH,  Ar,  6. ,  option,  R^ 

Here, there  are  five  items  to  be  input,  with  the  first 
two,  AH  and  Ar,  already  described.  The  third  item  is  6^ 
which  enters  (47) .  The  fourth  item  is  an  option.  If  the 
option  is  set  1  then  (48)  is  used  for  the  transform  of  the 
current  distribution.  If  the  option  is  o  then  (51)  is  used 
for  the  transform  of  the  current  distribution.  The  fifth 
item  is  conduction  loss  which  is  now  to  be  input  at  this 
point. 

Card  6-  heading  for  plots 

This  is  the  first  line  for  the  heading  for  several 
plot  frames.  Columns  1-70  may  be  used.  There  is  no 
indication  now  for  an  array  of  normal  modes. 


PSCGRAN  FO  CT  ( I NpL  T»OL  TPLT) 

CIMENSICN  PXCGIO  (3  ) 

CIMtNSICN  F< 12 jG ) ,FM(12u3) ,FP (1203 >,CAP  (12  GO) .CAM  (1200) ,VGM(120G> 
XFP(  1ZCQ  )  ,PM12  j  0)  ,KP(  12001  ,KN  (12u  31  «nU£CC)  ,  fX  (12CC)  , SERF (1200  )  , 
XSf.  Rf*(12U0l  ,V  NM  (  1 2  0  0 1  *  V  GP  (1  2 1  0  ) 

CiHENSICN  Hr!  flu  I  7  I,  PEA  LI  (7 1  ,H£AC2(7»  ,FEACi  (?) 

CIMENSION  FNC5G).VM50) 

CCHNCN  EL,  AL 1,  A  L  e  ,  5  *D  ,  Tl»G*S,ETA,EN,F  ,A  Y  ,  fl  »E  L  1  (40  I  ,FE  (40  ,  AA (4  0  I 
X .LMCCc ,CEL1, JOPT 
IzAC  *,H,Tl,C,G,EL,EN,ETA 
SiflC  *,EL9EGN,t LCEL.ELOPT 
-cflC  *,A«3ECIN,  ADEL,AOFT 
REAC  *,P9EC-IN,  PCEL.POFT 
FEAC  *,CFLH,DiST,CELl  ,  JCPT.FL 
PEflC  1JJ,HE40 
fi-AC  1 0  0  »H  EAC1 
PtflC  1 0  0  » H  £A02 
PEAC  1 CC ,H£AC3 
LN  C  C  E=  1 

IF  IPEAuJll  .EG.  "NCRf*  NCDE  “)  LM0DE=2 
N-Ffv 

EC  Rl  1  =  1, N 

ELI ( I ) =cl0EGN+ ( 1-1 )*EL  DEL 
flC(l)=A3EGIN4(I-ll*AOtL 
Pc(I)=P?EGIN+{T-l)*HOfcL 
N£L  =  (N+l)/2 

IF  (ELCPT  .EG.  0.)  GO  TC  h2 
nc  s3  I-NEL.N 

ELI <I)=EL1  (NEL  >  -  ( I -NEL  ) *ELOEL 
IF  (  flCPT  ,tC.  L . )  GO  TC  *4 
CC  U5  I  =  NE  L  ,  N 

AM  I)  =  AA  (NELI-  (  I-NEL>  *ACEL 
IF  (FCPT  .EO.  0.)  GO  1C  *6 
CC  U7  1=  NE  L  »  N 

Ft(II=PECNEL)~ ( I -NEL) * PCEL 
C  CMIMJE 

FlC  =  2.6*SQFr  (H*  (F  +  175U.M 
FH=2.8*  (H  4£ 75  . ) 

FCFL=1. 

FEECrfllM(FLC)  +1. 

NF=1NT  ( F  HI  )-INT(FLC)-l 
CC  43  1=1, NF 
F  (I  )=F3EG4  (1-1  >*FDEL 

IF  4 F  <  1  >  .LT.  FLC)  PRINT  *,  "FREQUENCIES  ICC  LCW" 

IF  (FtNF)  .GT.  FFI)  PkINT  *  ,"F  REQUENC I6S  TCC  HGH“ 

FRINT  60 
F  A  C  T  =  1  . 

IF  (ETA  .GT,  -2.  )  C-0  10  * 

E1A=-1. 

F AC  T  =  <2. /EM**  2 

PRINT  *,"  PI  CRATING  CASE” 

CONTINUE 

F-FIM  61, H  ,T1,  D,G,EL  ,EN  »  ETA  ,NF 

PRINT  *,"  DELTA  H  =  ",CELH,"  DISTANCE  =*‘,CI5T 

PAINT  *,  "  R  LCSS  IS  M  ,  RL 

F  T  I  NT  *,"  DELTA  =  *•  ,  DELI 

F  5  T  N  T  ei,(I,Eul(I)  ,1  =  1, N) 

PRIM  82,  (  I,  AA  (I  >  ,  1  =  1,  Nl 
PRINT  83,(I,PE(I),I=1,N) 

PRINT  feo 

TF  (ETA  .EC.  -1.)  ELL  =  ,87*  .4t  -6*  EN 
IF  (ETA  .EO.  1.)  ELL=  .23*.  4E-6/EN 
ELL  =  G . 

CAT B  FROGIC/dH SETH  ARES , 4H3724 , 10H J. HEINE  ERG/ 

CALL  FLTI33(PR0C-IC  ,20  0  •  ,12,  ,1  .  > 

P  =  0  . 

A  =  0  . 

A >=  w  . 

FI  =  2.lm59  2fcE*4 
RCMO. 

A  Y  T  =  «  5  *  A  Y  ♦  l(l.-ETA)+EN*U.tETA  I) 

LC=  4.*PI*1  •  E-7 


=  H/l 750  . 


CN  =  EF/ <2.oM75C.  ) 

L 11  =  1.-0MH/(CM**2-CMH  «*2> 

U22=Oli 

L12=CM/ (0M++2-0MH*«2) 

==SC5T (U11/U2?) 

L  =  1 

IF  <1  .EG.  2)  GO  TC  2 
S  =1 . 

GC  TC  3 
S  =  -i. 

CONTINUE 
A  L 1 =U2  2  + S*  S  '  11 2 
AL2=U22*C- S*L12 
IF  CJ  .GT  .  1)  GC  TO  b  3 

CAO=.5»SORT(12  2/U11)*ALOG(1.*<4.»SQRT  (L11«L22>  / 
X  CL12#*2-  (SORT  (U1 1*122  )  +1 • I ♦*  2  )  ) 

C AO=CAC/C 
CONTINUE 

IF  (I  .EC.  1)  GO  TC  51 
If  (J  .EC.  i)  GC  10  51 
I  F  <  L  .EQ.  1)  CAC  =  CAP< 1-1) 

IFCL  .EC.  2)  CAC=CAM(J-1> 

N  =  1 

CFL  =  .C2*CA  C 
C  ACF=CAO  +  D  Ei. 

C ACN=CAO-DEL 
C  AC C=C A0*0 
C  A C  C  =  A  3S  (CA0C1 
C  CO G  =  A 3S  <C  AO *G  ) 

I F  (  C  ACD  .GT.  o  5  0  •  )  GO  TC  35 
IFiCACG  .GT.  o  5  C  .  )  GO  TC  35 
F1CC=FT (CAO) 

FTCF=FT<CAOP) 
f  TC  F  =  FT (CACM) 

CA1=CCC-2.*0EL4  FTCC/l FTCP-FTCM) 

IF(A?E(CA1)  .GT.  1.E7)  GC  TO  35 

CAiC=CAl*0 

C  A 1 C  =  A  cS  (C  Ait) 

C  A 1 C  =  A  3S <CA1*G> 

IF (CA1D  .GT.  b 5 C  •  )  GO  TO  35 
IFUA1G  .GT.  65C.1  GO  TC  35 
FTC1=FT(CA1) 

IF  (A3S(  (CAl-C  AC1/CA0)  .LT.  .301)  GC  TC  10 
CA0=CA1 

N  =  M4l 

1 F  (  N  .GT.1G)  GO  TC  35 
GC  TC  5 

IF  (L  .EO.  2)  GO  TC  20 
CA=CA1 

IF  lArS(FTCl)  .CT.  1.)  GC  TO  35 

IF  <CA  .LT .  0. )  GC  TO  35 

FF  (  I )  =EF 

CAP  (I )  =  C A 

1=141 

L  =  2 

C-C  TC  3G 
C  A  =  C  A1 

IF  (A5SCFTC1)  .GT.  1.)  GC  TO  35 
IF  (CA  .Li*.  3.  )  GC  TO  35 
FM  J)  =  EF 
CAM  (J ) =C A 
0  =  J  4  1 

IF  ( J  .EQ.  I)  GC  TO  15 
IF  (J  .GT.  I)  GO  TC  31 
I  =  J-1 
K  =  K-1 
J  =  J-1 
K  =K  +1 

IF  (K  •  LE  •  NF)  GC  TO  50 
PRINT  60 
11=1-1 
Jl= J-l 
GC  TC  2*. 

PRINT  ♦.•'ITERATION  GOES  NOT  CONVERGE. F=  **  *  EF  »  * 
IF  (L  .EC.  2)  GO  TC  lb 
1=2 

GC  TC  2 


2t  CONTINUE 

PRIM  63,  (  FP  Cl  )  ,CAF(I  )  ,1=1*11,  It) ) 

PRINT  o4 , ( FM ( J )  , C  AN ( J  )  ,J=i,Jl,13> 

PRINT  60 
CC  22  J=1,J1 

IF  (J  .EG.  1)  VGN(«.)  =  5./FI*(CAN(2)-CfMll)/FCEL 
IF  (J  .EC.  Jl)  VGM  (J)  =5. /PI*  <CAM(  JD-CAMJ1-1  ) ) /F  CEL 

IF  (J  .NE.  i  .ANC.  J  .NE.  Jl)  VGM(JI= 

X  5  . /FT* (CAM'J*:  >-CAM(  J-i>> /PC  El*.  5 
22  CONTINUE 

PRIM  85,(FM(J)  ,VGMJ).J  =  1»J1»10) 

PRINT  60 
CC  21  1=1,11 

IF  (I  .tO.  1)  VGF(I)=t»./PI*(CAFC2)-C«F(l))/FGEL 
IF  (I  .EC.  ID  VGP  II)  =  5. /PI*  CCAPt  ID-CAF411-1  ))  /  F  CEL 

IF  (I  .NE.  1  .ANC.  I  .NE.  ID  VG»(I>= 

XE  •/ F I  *  (CAP  ( I  ♦!  )-CAF(I«l))/F0EL*,5 
21  CCNTINUE 

PRIM  57,  <F°  (I  )  ,VGF(I  J  ,1  =  1,11,  10  > 
j  C  =  Jl/2 
CAM" =CAM { J  C  ) 

CC  23  J=l,Ji 
C  AMN  =  CAMG- CAM(  J  ) 

C AMN=CIST*CAMN/PI 
fAMN=AMCCtCAMN,2.) 

IF  (CAMN  .GT.  1.)  CAMN  =CAMN-2 . 

IF  (CAMN  .IT.  -1.)  CAN  N  =  CA  MN ♦ 2  « 

VNN  (J)riPJ  ,  ♦  CA  MN 
23  CCNTTNUE 
PRIM  60 

PRIM  66,  <  F  M  ( J  )  , VNN( J  )  ,J  =  1,J1,1C) 

S  =  1  . 

CC  25  1=1,  11 
E  F  =  F  F  ( I ) 

C  A  -  C  AF  (  I ) 

C  =  2.*FI»FF*i,“f 
CN=EF  y (2.6  *1760  .  ) 

1)1=1. -CMH/C OM*  *  2- CMH*  *21 
L  2  2  -1)1 1 

112=CN/(0M**2-0MM*«2) 

==SCRT (U11/J22) 

A  L 1  =C2  2*  B+  S*-  UA  2 
A12=L22*B-S  U12 

Fl=  (T  (CA)*  i,  )♦  *  2  *  ( 1  ANN  (CA*  El)  -  CA*  EL*  I SECMC  A*  EL  1  1**2 ) 

F2=  (Rl (CA) *EXP (CA*G»-R2(CA»*c>F(-CA*G)> **2/4. 

X*  (CCTH(CA*TD-CA*T1*(CSCMCA*T1M  **21 
F2=  .25*  <Ri  (CA)  **2)  MEXP<2.*tA*C)-l.)-.2  5MR2(CA)**2> 

X*  (FXF(-2.*CA*G) -1 . )-Rl <CA)*R2(CA)  *CA*G 
Fm  =  AL1*(T(CA)**2»*  (tXP  (2.*8*CA*0)  -1 . ) -A  L  2*  (E  >  F  (-2  .*E*CA*D )  -1 . ) 
X-k.*'?**2*U22*CA*C*T(CA> 

G  EC  A  =  C£  1C  A  ) 

FF(I>=  0*'J0*  GECA  **2* (P1*F2**3*F4> ✓£  ./CA«*2*  .5 
FM  I ) =4 • *  P  F ( I) 

RFC  IJ  =  A  =  E(F°(I)  )/(  (l.-ETAI-Ml.  ■*ETAI*€N**2)*4. 

RF  <  J ) r FACT  *RP( I ) 

25  CCMIME 
5.  *  -  1 . 

lC  26  J= 1 , Jl 
EF=FN  (J) 
c  =  2  .*Fi*FF*i.F:e 
C  A  =  C  AN  (J) 

CN=EF/(2.o*1750.  ) 

111=1. -CMH /( CM **2- CMH *  *  2 ) 

122=111 

L  J.2-CN/  (OM*4  2-OMP*  *21 
C-ECR7 (U11/U22 ) 

AL1=U22*3*S  Cl  2 
Al?=L22*?-E*ll2 

<T  (CA)  M  .  V  *2*  (TANN  (C  A*  E  l )  -  CA*  E  l  *  (St  CMC  A*  EL)  »  «*2) 

F2  =  (Rl (CA) *FXP (CA*C)-h2(CA)*EXF(-CA*C >)  **2/4. 
x*  (CCTMCA*T1  >-CA*l  1*  (CSCMCA*T1»I  **21 
F*=  .25*  (Fl  (CA)**2)*(EXPC2.*CA*C)-l.l-.2  5*  (R2(CA)**2» 

X*(f >P(-Z.*CA*G) - 1 . >-R 1 (CA) *F2 (CA> *CA*C 
FM-ALl*(T(CA)t*2)*(tXF(2,*8*CA*D)-l.)-AL2*  IE>  F  (-2  .*S*CA*0) -1.  ) 
Y“‘-.*E**2*U22*CA*C*T(CA) 

G : C  A  =C  c  (  CA  ) 

FF(J)=  0*U0*  GECA  **2* CP1  +  F2+P3+FR) /6./CA**2  *.5 
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FN  (  J  >  =  APS  (  PM  (J  )  >  /  (  <1 .  -  ETA)  -Ml  .♦ETA)"EN""2)"4. 

RN  (  w)  =  FACT  "RN  <  J ) 
cl  CONTINUE 
PRIM  60 

F  R  I  NT  71  ,(  FP  (I  ),  Kp  (1)  ,  1*1,  II,  1  C) 

C  PRIM  72,  (  FM(J)  ,FM  (J)  ,  J  =  l,  Jl,5) 

IF  (II  •  NZ  •  -1)  GC  TO  90 
CO  64  1-1,11 
P.M  RI(I)=fiP(I)*RH(H 
PRINT  60 

PRIM  75,  (F°(I)  ,  FT  (I)  ,1  =  1,11,10) 

FF1=FF (1) 

F  FL  =FF ( II ) 

CALL  FTRAN  RT,  PX  ,1 1  ,FP  1,  FPL) 

PRINT  60 
N  =  NN 
M=I1-1 

PRINT  73  ,  (  FP  (I  )  ,FX  (I)  ,  1  =  1, M  ,20) 

CC  54  1=1, N 
XL=2.*FI*FP(I)"ELL 
Xt=XL"l.E6 
L  =  I 

IF  (I  « EC*  1)  l-l 

SERF  (I)  =  2C.«ALCGiQ((4."RP(l)/RG)/ 

X  (  (1.MR7  (I)  +  RL)/FGM*2 

X  +  (  (  PX  (I )  4-XL  )/30)**  c  I  ) 

SERF(I)=SER°  (I)  -76  .4*  1)  ELH*D  I  ST  / 

X(2.*PI*1.E6*  (FP(IM)-FP(L-l))/  (CAP  (I  #11  -  CAP  (L  -1)  )  ) 

SSRMI)*  2  0.*ALCC-1C(  (  R  .*RM  (I)/RG)  / 

X  (  (l.*(FT  (II+RD/RG)  **2 
X  +  (  (PX(I)*XL  )/RG>**<)  ) 

$€«?MI)=SEfM(I  )  -  76  .4"OELH*DIST/ 

X(2.  "FIM.E6-  (FM  U4l)-KM(l-l)  )/  (CAM(H  1)  -CAP  (l  -1)  )  ) 

£4  CONTINUE 
PRINT  60 

PRINT  77, ( Fc  (I )  ,SEFP( 1 ) ,1=1, N  ,10) 

PRINT  60 

PRINT  76  ,  ( FP  (I )  ,SERM( 1  )  ,1  =  1  ,M  ,10) 

XMN  =  2500. 

C  X  =  5  C  » 

XNIN=25GC. 

XNIN=240r . 

F  6E  GC= • 0 1*  Fa£G 
X  (•  I  N=  A  INT  (  FI  EGO)  "100. 

)>  =  AINT (.0 1  +  FHI) -AINT ( .01* FLC)  ♦!. 

C>=  ICC. 

YNIN=C. 

CY=200. 

CY=1CC0C . 

CC  96  Jsl.Jl 

IF  (CAM  ( J )  .  GT  .  1CC00  0  .)  C  A  M  ( J  )  =  1  3  0  0  C  C  . 

IF  (VGM  ( J )  .GT.  1CCC.)  VGM ( J )  =  1 0  0  J  • 

56  CONTINUE 

CC  57  1  =  1,  II 

IF(CAf(I)  .GT.  1GJCU0.)  CAP (I) =1C 303 0  . 

IF  C  W  G  P ( I )  .GT.  1CC0.)  VGP ( I )  =  100  3  . 

97  CONTINUE 

CALL  flCTCl.5,0.,-3) 

CALI  $YM°0L(.3,5.6,.l  ,HEAD,C,7G) 

CALL  SYN=>OL<  .-  ,5.6  ,.l  ,  HE  ADI  ,0 ,70  ) 

CALL  S  YM  PQ  L  (  .5 ,  5  •  h  , . 1 ,HEA02,U, 70) 

CALL  SYM50L ( .5 ,5.2  ,.l  ,HEAD3 ,3,70 

CALL  AXIS(  C.  •  ,  2  1  H*A  V  E  NUMBER  (♦  )  (1/M)  ,21,1(.,9C»,YMIN,0Y,1C.  ) 

CALL  AXIS  (C.,  0.  ,1£HF^ECLENCY  (MHZ  I  ,  -  15  ,  >  >  ,  C  . ,  XM  ]  N  ,CX  ,  10  .  ) 

CALL  LINE(FM,CAM,jl,l,0,l,XMIN,OX,YMIN,CY,.06  ) 

CY=i00. 

CALL  FLCT(  16. ,0., -3) 

CALL  A X IS (  0.  ,J  .  , 

X  25HGRCUP  DELAY /CM  (♦)  ( NSEC  )  ,  25 , 10 . , 9C  .  , YMIN , DY , 1 0  .  ) 

CALL  AXIS  (  i  .,  C.  ,1  EHFR  EQUENCY  (MHZ), -15, XX  ,  b  . ,  XM  I N  ,0  X  ,  10  .  ) 

CALL  LINE(FM,VGM,J1,1,0,1,XMIN  ,DX  ,  YMI  N  ,  t  Y  ,  .  08  ) 

CY=  ICCGu  . 

CALL  f  LCT  ( 16  • ,  0  •  »-  3) 

CALL  AXIS( C. ,0 . ,21FWAVE  NUMBER  (-)  ( 1 /M  )  ,  2 1  , 1 C .  ,  5 C  .  .  Y MIN, 0 Y , 10 .  ) 

CALL  AXIS  (-.,C.  ,15HF*ECLENCY  (MHZ), -15, XX  ,  C  . ,  X  K  I N  ,C  X  ,  10  .  ) 

CALL  LINE( FP, CAP, 11,1 , 0,1, XMIN  ,DX , YMIN, CY  ,.06  ) 
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C  1  =  4  L  >.  I 

can  plc7( ie.,c 
ran  axis<g.  ,c. , 

X  ?  5  HGRCUP  DtLAY/CM  <- )  ( K St  C  I  , 25 , 10  .  ,  «G  . ,  YH1N ,  OY ,  1  fl  .  ) 

CALL  AXIS  (,  0.  ,  INFREQUENCY  (MHZ),-15,XX  ,  C  .  ,  XM  IN  ,CX  ,  10  .  ) 

CALL  LZNt'fF'tVGF  ,11,1, 0 ,1 ,  XMIN  ,OX  ,  YMI F  ,  CY  ,.0E  ) 

YCIKs-160. 

0=36. 

CALL  FLCT ( lo • , 0 • , • 2 ) 

CALL  AXIS< G. ,0. ,<2HNORMAl  C1SFERSI0N  (♦>  ,22 ♦ 1 0 .  ,  5 C  .  . YMIN , OY , 10  . ) 

CALL  AXIS  <  J  G.  ,15HF*ECU£NCY  CMHZ),-15,>X  ,0  .»XH1N,CX,1G.) 

j  j  =  n 

5  ?  CCNTTNUE 

LC  59  J=1.J1 

IF  UJJ  +  J)  ,  GT .  Jl)  GO  TO  57 
FN  <  J)  =  FM  (  J  J  +  J) 

VM  <L)=VNM< ^J+J) 

IF  (J  .EG.  1)  GO  1C  59 
IF  (IM(J)  .LE.  VM(j-l))  GO  TO  59 
GC  TO  57 
5C  CCMINL-I 

c  j  '  -  v  ■ 

CALL  LINE  (FN.VMtJJ,l,o»l,XMIN,OX,YMIhf  Cl,.Cfi) 

*...  =  JJ-t  J2 

IF  <  (JJ*1  )  .L  £.  Jl)  GC  TO  50 
YM I N  =  0  . 

0=20. 

CC  27  1=1,11 

IF  (KF(I)  .GT.  3  C  Q  . )  HrP(II=30Q. 

IF  (  SP  1 1 )  ,  GT  •  2GC0.)  RP(I)=23oC. 

27  CONTINUE 

CALL  FL  CT  (  15  •»  u  .,-2) 

C  a  L  L  AXIS  .,0  .  ,27HRflO,  RES. .MINUS  tiAVE  < CHM  S)  ,  2  7 , 1 0  .  ,  9u  .  , 

x  ymin,cy,i:.> 

CALL  AXIS  ( j  . ,  0  .  .  1  EHFR  ECLcNCY  IMHZI,-15,XX  .0 . , XM  IN  ,CX , IQ  .  > 

CALL  L INl ( FP»RP,I1 ,1, J «1«XMIN, CX,YMlN,OY».0*) 

CC  28  J=1,J1 

IF  (RM(J)  .GT.  3Co  •)  RM{J)=30«. 

IF  (FM(j)  .GT.  2CCC.)  RMJ)=2:3C. 

IF  (KT(J)  .GT.  3  CO  .)  RT(J)  =30-  • 

IF  (R7(J)  .GT.  2C0C.)  RT(J)=20CQ. 
if-  C  CM  TINGE 

CALL  FL0T<17.,  C  .,-3) 

CALL  AXIS  (j.,G.  ,2  7HR  A  0*  RES. ,  PLUS  RAVE  (CHM I) , 2  7 ,10 . ,50  . , 

X  Y  M I N  ,  L  Y  ,  10.) 

CALL  AXIS  CJ  .,  0  . ,  1EHFR  EQUENCY  (MHZ1,-J5,>X  ,  0  ,  ,  XH  J  N  ,CX  ,  10  .  ) 

CALL  LINF(FM,PM,Jl.l,0,l,XMIN,CX,YMlN,OY,.ufc) 

CALL  FLCT(  17. , 0. ,- i) 

CALL  AXIS(  0.  ,2  .  ,  22HR4L  .  RES.  TOTAL  (CM*S  )  ,22, 10.  ,50., YMIN,  OY  ,10  .  ) 
CAL*.  AXIS  (1  .,  0  ,  ,1  5HFy  tQLENCY  (MHZ),-15,>>  ,0  ,»XM  IN, LX, 10.  ) 

CALL  lINE(FP,PT,I1,1, U ,l,XMIN,CX,YMlN,DY,.Go) 

1 1  =  M 
J  1  =  M 

YMha-lOCO  L. 

YMIN=-250. 

o=2C:o. 

0  =  5  j. 

CC  ;  2  1  =  1,11 

IF  (PX(I)  ,LT.  -1CCL0.)  FX(I)  =  -10J0C  . 

IF  (FX(I)  .LT.  -25C.)  PX(I)=-25u. 

IF  ( F  X ( I )  .GT.  1CCC0.)  PX(I)  =  1C00  0. 

IF  (FX  ( I )  .GT.  25  C « )  FX(I)=25C. 

52  CONTINUE 

CALL  r  L  CT  <  17.  ,C  .  ,-  2) 

CALL  AXIS(  C.  ,0.  ,  22FRAL  .  REAC  TCTAL  (  C  FMS  )  ,  2  2 , 1 0  .  ,  50  .  ,  YMIN  ,  OY  ,  10  .  ) 
CALL  AXIS  (0  .,  0  .  ,15HF(i  ECUENCY  CMHZ),-15,>X  ,C  .,  XM  IN  ,CX,  19  .  > 

CALL  L INF(  FP,PX , 11 ,1, 0 ,1,XMIN,CX, YMIN  ,OY  ,.06) 

CC  c?  1=1,11 

IF  (SERP(l)  .LT.  -80.  )  SERP  <  1>  = -8  j  . 

IF  (SERM(I)  .IT.  -80.  )  SERM(I)=-80. 

53  CCMTNLE 
YMIN=-Cj. 

D  V  =  1  G  . 

call’flctc ie.,l .  ,-2> 

CALL  SYMBOL!  .  J  ,5.0  ,  .1 ,  HtAO  ,  0 , 7  0) 

CALL  SYM°0L (  .5 , S  .  6  , . 1 , HEAD  1 , 0 , 70 ) 

CALL  SYMBOL(.c  ,  9  . 4  , .  1 ,  HEAD  2  ,  u  ,  7u  ) 
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CALL  SYMBOL ( .5  ,  S  .2  ,  .1 ,  HEAD3 . 0  ,  70) 

CHL  AXIS  (  C. , - .  ,  261 — IN  S.  LOSS, MINUS  WAVE  (C9>,26,  2.,SC.,YMIN, 

x  nv,io.) 

CALL  AXIS  (J.,C.  , 15HF*£CUENCY  (MHZ), -15,  XX  ,0  XH  IN  ,CX  ,  10  .  ) 

TALL  LINF(  FP,SLRF,I1,1  ,(J,1  ,XMIK,9X,YMlN,CY,.0f) 

CALL  FL  CT  (  16  . ,  C  .  ,-3) 

CALL  SYMEOLl.5,5.6,.l,hEAD,0,7ti> 

CALL  SYMBOL  <  .5 ,9  .fc,.l  ,HEAQ1,  (,  70) 

CALL  SYM80  L  (  .5 , 9  .  *,  , .  1 ,  HE  AO  2 , 0  ,  70  ) 

CALL  SYM?OL(.5,5.2,.i,H6A0  3,0,70l 

CALL  AXIS<  C.  ,9  .  ,26H-INS.  LCSS,  **LUS  WAVE  ICE)  ,26,  t.,90.,YMIN, 

X  0 Y, 1 _ , ) 

CALL  AXIS  (u.,  0.  ,15HFREQUENCY  (MHZ), -15, XX  ,  C  .,  XM  ]  N  «C  X  «  10  <  ) 

CALL  LINECFM,S£RM,11, 1  ,  0 ,1 , XMIN, OX , YM  IN  , C Y  ,  .  C £ ) 

CALL  ENDPLT 
STOP 

PRINT  76 
STOP 

FCRMAT (1H1 ) 

FCPMAT(5X,"  H=  **  ,  E 15  .  7  /  5X  ,  ”  Tl=  **  , El 5  .7 /5X  ,**  C=  **,El5.7/ 

X  5  X  ,  **  G=  **»  ci  5.  7  /  5X  .**  L  =  **,£15.7/ 

X  5  X  ,  “  N  =  **,  Ei  5.  7  /  5  X  ,**t  T  A  =  '*,  El5. 7////B  X  ,**  N  C  •  CF  F*S  AFE  **,I5) 

FCRMAT (//ICC 5E 15*5/)) 

FCRMAT  (  /  //  51 X,  **  S  =  1*V  /  <1JX»**F  =  E 15 . 1  ,1 0  >  ,"K  (-)  =  **,El5.7/>> 

FCRMAT  <///5jX,  **S  =  -1*V/  <1QX,**F  =  ** ,  E15  .  / ,  1  L  X  ,**K  (♦)  =  ".E15.7/D 

FCRMAT  (///  (1CX,**F  =  “,  E  IS  .  7 , 10X  »“P  (-1  =  **,E1E.7/)) 

FCFMA7C///  (1CX,**F=  **, 615.7, 10X.**P  ( ♦  )  =  **,E15.7/)) 

FCRMAT  ( /  /  /  5  X  ,**  L  1=  "  ,E  1  5 . 7/ 5X  ,**  A=  ** ,  £  1  5.  7  /  5  X  ,"  f =  ",£15.7/ 
X5X»“IC=  **,  Ei5.7/5X,**N  =  "  ,  El  5 .  ?  /5X  ,  **E  T  A  =  -,E15.7) 

FCRMAT  (///  (luX,*‘ F=  **,cl5«7»lJX,**RAO.  Ft  S  .  (-)  =  “,615.7/)) 

FCFMAT  (/✓/  ( 1  OX  ,*’F  =  “,tl5.7 ,10X  ,**RAO.  FES.  (♦)  =  “,El5.7/)) 

FCRMAT(  //  ( IQX  ,  **F  =  **,tl5.7,lOX  ,**RAD.  FEAC.  TC  TAL=  “,E15 .7/)) 
FCRMAT  (/✓/  (10X  ,  *'F  =  ** ,  1 15 . 7 , 1 C  X  ,**R  A  0  .  FES.  TCT  AL  =  ",El5.7/)> 
FCRMAT  ("1  A  K  RCCT  tXISTS  FCF  ONE  WAVE  CNL Y “) 

FCRMAT  (//(  lgX,**F=  * ,  E  1  5. 7 , 1C  X  ,  "INS .  LCSS  (->  =  ",El5.7/>> 

FCRMAT  (//(  i>  X,  *' F  =  ,E  1  5 . 7 , 10  X  ,  “IN  S  .  LCSS  (♦>  =  **,615.7/)) 

FCRMAT  (//(  Ij  X,**L1  ( “ ,  I A  ,“)=  **,E15.  7/)  ) 

FCRMAT  (//(  lgX,"  A  ( *'  *  I A  ,*’)=  **,615.7/)) 

FCRMAT  (//(  UX,**  F  (",  IA  ,**)  =  **,£15.7/)) 

FCRMAT  (/(1CX,»  F=  ** ,c  1 5. 7 , 10 X  ,  **GR OUP  tf  LAY  («>  s  “,615.7/)) 
FCFMAT  (**  F=  ** ,  €  1 5  •  7  ,  '*  NORM  CISPERSTON  (♦)  =  **,E15.7) 

FCRMAT  ( **  f*  **,615.7,**  GRCUP  DELAY  (-)  s  **,615.7) 

FCRMAT (7Alt) 

EMC 


FLNCTION  COTH(CA) 

CCMMCN  EL, AL1,AL2,R,0,T1,G,S,ETA,EN,F,AY,A 
CCTM  =  1./TANH  (CA  ) 

R ETLFK 
£N0 


FINCTICN  SECH(CA) 

CCMMCK  EL,Aul,AL2,F,D,Tl,G,S,tTA,EN,F,AY,A 

IF  (CA  *  •  Lt  •  7w(3.)  StUH  =  l./COSH  ICA  > 

•*?TLRN’ 

END 


FLNCTION  CSCH(CA) 

CCMMCN  EL , AL 1 , AL  2 ,E*0*T1*G,S»STA,EN,F  ,  A  Y  ,  A 
CSCMsC. 

IF  (CA  •  L£ .  7*0.)  CSGH=1./SINMCA) 

RETURN 

END 
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F INC7  I  ON  hTUCA) 

CCP*CN  £L»ALl»ALZ»?»n»7lfGfSt£.1A*cN*FtAV,A 
F  1 7  =EXP<-9*CA*C)* (Ri(CA)*EXP(CA*C l-R*(CA)«EXF  l-CA*G)> 

X«SM1MCSCMCAM1> )**c 
F  7  7  2  = 

X-E»C*EXP(-?*CA*C>*  <<COTMCA*Tl  14-1  ,)*"U  CA)«E>  F(CA*0» 

X4  «CCTMCA*Tl)-l.)*F2tLAI*EXP<-CA*G)  ) 

F  77  2  = 

X4a.i:7*S*0*EXP{-2.»5*CA*C)*  (  CCC7H(CA*71)  4 1 .  I*  ( 1  .-A  12 1 
X*EXFICA*G)  -<COTMte*Tl  (i  .4-AL2MEXF  l-CA*  O  » 

F  7  7  5  = 

X  ♦  S  *  EL*  (  A  L 1  ♦  A 12  )  *  (  SECH  C  CA*EL  )  •*  2)/  CALi-TAht-  <CA«EL)  1**2 
F  7  7  t  = 

X  ( <CC 7M CA  #Tl)  -1.1  <71  . -All  >  *  EX  F  C-CA*C  >-  ICC7M  CAM  1)41.1 
*M1 .4Al17*£XP(CA«G>> 

F77A=F7T5*  FTT6 
F71=.E*(FT74FTT24F7734F7Th) 

SF7LFh 

EKD 


FINC7I0N  T (CA) 

CCMCN  c.i»All,Al2*c*D»Tl»G*S»t7A,ENtFjAY.i 
I  =CAL2  4TANH<CA*tL  )  )/IALl-TANH  (CMEL  M 
RE7  LFN 
END 


f  i  >,r,CN  =  ;5»  UCOTMC  A»71  )-l^»*R27(Ci  A*G>  «EXF  «-e*CA*D> 

X-lCC7MCA*Ti)  +  i.)*filU;A)*EXPCC  A*G  )*E>M-?*CA*CI> 

PET  LF^ 

END 


FLNC7I0N  GE(CA) 
RETUFN 

c  k  n 


AAUU  ) 


FI. NC 7 T ON  SINC(CA) 

Pl^JiA159£6  5ilAl  2  f  E»°* 11  »G»S»E7At  EN,F  ,AY  ,A,EL  :<40  ),FE(40)  ,AA  (43  > 
SlNf=<SIN(PT*r.  A)  »  /  ip  t  ♦  r  a  * 


FLKCTICN  GAY(CA) 

CCMFLEX  DC  ,D1,C1,F1,A1  ,A2,A3,A<,,C,CS 

xaSccE  !flilj5p^?’£*d,71,g,s’eta*£n,p’av,#*EL3<<4Q  >*FE<4ii)  »AA<*»0  > 

P 1=  2  •  1  **1FS  2c  Ew 

MEN 

MN 

IF  ( JCPT  ,  EG ,  i)  M=i 
£  =C NFLX ( c . ,:.) 

CC= .5*0EL1*CPPI X  (1  .,-1  . ) 

£C  1  1  =  1, M 
MAA (I) 

El  =  ELKI) 

Fl  =  FF  CI> 

C  1=  •  £  *  H  /  C  C 
ClR=FcAL (01) 

C  1 1* AIMAG  t Cl) 

£1=  .5*WMC  A  +  l./CC) 

F1=.5*WMCA-1,/CC) 

A 1  =  Cl-  (CCOS(Cl)+CCOS  (Fi)  >  +  Fi*CSIN(C  1)  ♦CSIMF1)*C1 
A2  =  -ClMCC0S(Cl) +CCOS  (Fl>  >+Fi*CSIN(Cl) ♦  tS!MFl)*Cl 
A  2=  CiMCSlN(Fl  ) -C ' IN  (Cl)  )  +  F1*  CCOS  (Cl  )-Cl*CCC5(Fl) 

A*-=  Cl*  (CS  IN  (Fl)-CSIN  (CD)-Fl-CCOS  ( C  1  > ♦  C  1*  CCC  S  <  F 1 ) 

CS  =  CFPLX  (C  CS  (CA*F1*I)  ,  -S  IN  (C  A  *  PI*  1 1)  *  El  A**  I 

IF  (  JC^T  <cC.  1)  CS  =  ( 1  •  -ETA**  N*  CMPLX  (CCS(CA*Pl*EN)  »SIMCA*Pl*EN)  )) 
X  <i.-ETA*CF,FLX  (CCS(CA*P1>  ,SIN  (CA*P1  )  )) 

C  =C  ♦SCPT  (El)  *CS* 

X(Al*CEXP(01MA2*CcXP(-Dl)  +  CHPL>(0  .  ,  - 1 . )  ♦  (  A  2  *C  EXP  (  C 1 )  ♦  A4*Cc  XP  (-0 1 

X  (W**2MCA*  2«-C)'FlMl.  ,1.)  *2./  CC**2)  * 

X  (CHFLX  (  COS  (Cl  I  )*S  1NH  (Olfi)  ,SIN  (D1 I)*CCSMC1F>  )/Dl  ) ) 

CCNIINOE 
G  AY  =C A  PS  (C  ) 

RETLRN 

ENC 


FINCTIOK  Rl(CA) 

CCHFCh  EL,  AL1,  AL2,  S  ,0  ,  U  ,G  ,S  ,E  1A  ,  EN  ,  F  .«  Y  ,  A 
U  =  <1  .-AL2)  *cXP(-S*CA*0>  ♦  II.  *AL  1)*  T  (CA)*  EXP(e*CA*0) 
FETLRN 
ENO 


FINCIICN  R  2 ( CA ) 

CCHPCN  £L,AL1*AL2,P,0,T1,G,S,£IA,EN,F,AY,A 

R  c  =(1.  ♦  U2>*EXF(-“*CA*0>  +  (1.-  AL1)«T  (CA)*E  >P<8*CA*C> 

RETL5N 

END 
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Reproduced  from 
betl  available  copy. 


S13PCLTINE  HTRAN  <R  ,X,  f  ,F3£G,F£NQ> 

CiHEKSION  R(3),XC3) 

PI=3.1i*l5926535E 
FC  EL=  (FENO-FEEG)  /  (N-ll 
F  =  F  =  EG4.5*FQEL 
1NC=NPD <N, Z) 

M-MTNC-1 
FH1=N-1 
MH2  =  M-2 
CC  23  1=1, NM1 

x  m  = : . 

IF  (I  .EG.  1)  RX=J3.*R  <1>*6.*R  (2)*R<2»)/£. 

IF  (I  .EG.  NNi )  RX  =  «-H  (N-2)+6.«R<NHlM3.«R  (R)>/8. 
IF  (I  .EG.  1  .OR.  I  .tQ.  NH1)  GO  TO  2C 
RX=  <»KCI-l)^S.*R(I)+9.*RII+i)-RCI+2ll/16. 

20  CCNTINUE 
F  1=  F  PPG 

CC  2b  IF=1,NIM2,2 

X  <1  >=X  (I)  +  4.*(R(IF41)-RX)/U  FI*FDEl>«*2-F**2) 

X  +  2.  *  ( R  (IP  I  -  RX)  /  f  FI  **2-F**2> 

F1  =  FI+2.»FCFL 
Zb  CCNTINUE 
F  EN  =  F ENG 

I F  (  INC  .EQ.  0  )  FE  MFt  NC-FDEL 
>  (I)  =  X(I)+  {\tNI)-R>)/(FEN**2“F**2) 

X  -  <R  <i  J  -RXI/I  F8EG**2-F*»2> 

X  ( I )  =  FOE  L/  3.  *X  ( I ) 

IF (  INC  . EQ  .1 )  GC  7  C  30 

X  CIl*xm+  .5*FDEL*  UR(NII-RX)/  IFE  N**  2  -F*  *  2  ) 

X  ♦  IR(N)-RXI/<  FENC**2-F#*2>  I 

2  f  X ( I ) =2 • /PI *F*Y  (  I )  ♦  RX/  H*ALGG 

X  <  < l.-F/FENO)/  (l.  +  F/FtNCIMF+F  EcG  >  /  (F -F  E£G )) 

F  =  F  4FCEL 
33  CCNTINUE 
NN2-N-2 

Xl=  <  1 E  .  * X  «  1)  -10  .  *X  (2)  +  3.»X  ( 3  >  )/6. 

X2  =  C2.*X(l)+6.*X(2l-X(3n/e. 

CC  31  1=3, NM2 

X7  =  <-X  (I-2l  +  9.*XII-l»*9.*X(l)-X(I  +  lll/lt. 

X (I-2)=X1 

XI  =  X2 
X  2  =  XT 

31  CCNTIME 

X  tN)=<l5.*X<NMl)  -1C.*X  (NH2I  ♦  ?  .*X  IN-3  I  )✓  c  . 

X  <N-1)  =  (?,  *X<NMl>46.*MNM2)-.XCN-3»  )✓«, 

X  <N-2 )  =X2 
X ( N - 3 ) =X  1 
RLTLFN 
FNE 


Figure  1  -  Geometry  of  the  System 


Figure  2  -  Hyperbolic  Current  Distribution 
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Figure^  3  -  Dispersion  Relations  -  Generalized  Forward  Volume  Waves 
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Figure  4  -  Dispersion  Relations  -  Generalized  Backward  volume  Waves 
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Figure  8  -  Dispersion  Relation  -  Forward  Volume  Waves 
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FREQUENCY  l  MHZ  )  l  XI O' 


FREQUENCY  (MHZHXIO1) 
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Figure  \A  -  Group  Delay  -  Backward  Volume  Waves 
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Figure  17  -  Insertion  toss  -  Backward  Volume  Waves 
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